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THERMAL ELASTOPLASTIC ANALYSIS OF 


NON-METALLIC THERMAL PROTECTION SYSTEMS 
T. J. Chung and C. Yagawa 
ABSTRACT 

This report presents an incremental theory and numerical procedure to 
analyze a three-dimensional thermoclastoplastic structure subjected to high 
temperature, surface heat flux, and volume heat supply as well as mechanical 
loadings. Heat conduction equations and equilibrium equations are derived 
by assuming a specific form of incremental free energy, entropy, stresses 
and heat flux together with the first and second laws of thermodynamics, 
von Mises yield criteria and Prandtl-Reuss flow rule. The finite element 
discretization using the linear isotropic three-dimensional element for 
the space domain and a difference operator corresponding to a linear variation 
of temperature within a small time increment for the time domain lead to 
systematic solutions of temperature distribution and displacement and stress 
fields. Various boundary conditions such as insulated surfaces and convec- 
tion through uninsulated surface can be easily treated. To demonstrate 

effectiveness of the present formulation a number of example problems are 
presented . 


v. 
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1. INTRODUCTION 

The mechanics of thermoelastoplastic solids has attracted the attention of 
many investigators in recent years. The development of the theories of thermo- 
elastoplasticity began with an attempt to consider the plastic deformation 
as a thermodynamic state variable and with the controversial treatment 
of finite deformations. Crucial difficulties lie in the proper choice of 
free energy functional and methods of numerical computation. 

Biot [1] and Coleman and Gurtin [2] used the concept of state variable 
in dealing with thermodynamics of viscoelastic materials. Perzyna and 
Wojno [3"], Kratochvil and Dillon [4], subsequently, employed the similar 
concept by decomposing the total strains into elastic and plastic strain 
components. The recent work by Valanis [51 uses also the concept of hidden 
variables, but introducing no yield surfaces. Finite elastoplastic deforma- 
tions with thermodynamic considerations were also discussed by Green and 
Naghdi f 7 “] . Schapery [8] studied a thermodynamic constitutive theory with 
history effects represented by single integrals. Recently, Oden and Bhandari 
[9] proposed a general functional theory of thermodynamics of viscoelasto- 
plastic solids and showed that various theories of viscoelasticity, rate 
independent viscoplasticity and plasticity can be obtained as special cases 
by imposing suitable constraints on the material parameters. 

All of these theories mentioned above, however, lead to considerable 
mathematical and computational difficulties in dealing with problems of 
geometrical complexity. The present paper is an attempt to establish an 
alternate approach. We derive incremental governing equations for three- 
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dimensional thermoelastoplastic solids and provide convenient solution 
techniques by finite elements in space domain and the difference operator 
in time domain. Although an extension to materials with memory presents 
no special difficulty, the present study is not intended for viscous be- 
havior. In this study, the Lagrangian coordinates are used disregarding 
the finite strains, but large deformations are included in the formulation. 

Finally, example problems are presented to demonstrate the merits of 
the present formulation. Effects of thermomechanical coupling, thermoplastic 
deformations and stress, and developments of plastic regions as a function 
of time are shown and evaluated. Various features not included in this study 
such as temperature dependent material properties or dynamically coupled 

inertia effects may be treated with minor modifications. 

\ 

. 2. THERMOMECHANICAL PRELIMINARIES 

It is assumed that the behavior of the body under thermomechanical 
environments obeys the laws of conservation of mass, balance of linear 
and angular momentum, conservation of energy, and the Clausius-Duhem 
inequality, 

Po dv o ~ 

'o 

ct ’j 3 + pF J • - pin = 0 , a 1 3 = o 11 
pe = CT 1 J y u + q J -!- ph 

I 

D 4-^ q'e,! £ 0 


(lb) 

(lc) 

(l d) 



I P d v 

J v 


(la) 
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D = p9.7] - q|^ - ph (le) 

Here p is the mass density with the subscript o indicating undeformed 
configuration. criJ is the second Piola-Kirchhof f stress tensor; superposed 
dots represent time rates; strokes and commas are covariant and ordinary 
differentiations; FJ and uJ are the body forces and displacements; Yl is 
the strain tensor; e, h and 7] are the internal energy, heat supply and entropy 
per unit mass; qJ is the heat flux per unit area; 0 is the absolute temperature; 
D is the internal dissipation. It is understood that for small strains 
p 0 = p and for rectangular cartesian coordinates covariant and ordinary 
differentiations are the same. 

The free energy cp is defined as 

cp = r - eii (2) 

which leads to 

Pep = CT 1J y u - D - pT|0 (3) 

3. T HE RMOE LA S TOP LA S 1 1 C BEHAVIOR 

As discussed in Section 1, our approach is to avoid the functional 
theory or the state variable concept. Instead of expressing the free 
energy as functionals of all the histories of strains and temperature and as 
a function of temperature gradient, we postulate an existence of (D in incre- 
mental quantity so that for a small time interval At 

cp(At) = (At), Y <J>(At), Q (At) } (4a) 

Similarly, 

tf 1J (At) = (At), Yj ^ (At) , 0 (At) } 


(4b) 



4 


q 1 (At) = QCvia 5 (At), Yi j ^ (At) , e(At)} (4e) 

/N 

11 (At) = H(yi (At) , Y <J>(At), 0 (At) } (4d) 

Here it is assumed that the total strain is the sum of the elastic and 
plastic components with («) and (p) representing, respectively, ’’elastic" 
and "plastic". 

The implication of (4a) through (4d) is that the free energy, stresses, 
heat flux, and entropy are functions of elastic strain, plastic strain and 
temperature only within the small time interval. It is then a simple matter 
to derive the governing heat conduction equations and equations of equilibrium 
in "incremental form". All histories may be carried over from one increment 
to another through the step-by-step numerical time integration. 

Although the free energy can be shown in various forms depending on the 
material properties or the purpose of analysis, the present study is limited 
to a particular case including only quadratic terms as follows: 


pep (At) = ^ilW Y ( ( »)(At)Y'||, t H- ^ ,J 'cV 1 5>(At) Y ^W 

-B ,J T(At) y (e>(At) - llUT(At) v <y>(At) 

c T g ( A t) 

2T„ 


(5) 


7C 

in which EiJM and EUkl are tensors of elastic and plastic moduli; B'J and 
A- 

B 1J are tensors of thermoelastic and thermoplastic moduli; c is the specific 

heat; T 0 and T are the reference temperature and temperature change related 

* * 

by 9 - T 0 + T. Here EHki and B'J are the functions of current state of 
stress in the inelastic range. It should be noted that such inelastic material 
properties cannot be admitted had the form of free energy been expressed in 
entire history domain rather than in a small time interval. These arrays of 
plastic moduli remain constant only during the small time interval and must 
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be updated as the state of stress changes. 

Likewise, the expression (3) in a small time increment is given by 

pco(At) = ct 1 3 (At) Yj j (At) - D(At) - pT)(At)9(At) 

= ffU (At) (y^ (At) + yi ^ (At) ) - D(At) - pTl(At)T(At) (6) 
and the time rate of (5) becomes 

pcp(At) = K Uk£ Y ( v2 ) (At); ( k « ) (At) + EU^ Y y>(At)y( k J>(At) 

- b’ J T(At) Y ( 1 « ) (At) - B 1 lT(At)y ( 1 j ) (At) 

- h JTCAtly^^ (At) - BUT(At)y< j’ - - — ^ (7) 

In view of (6) and (7) and dropping (At) for simpler notation, we obtain 

(Ei - BUT - ctU)^^ + (-BU y (j) - bU y (J) 


c T 


+ P T)) T + D + E 1 - a 1 J Y l j } = 0 < 8 > 


For all arbitrary values of Y <®) and T the following relationships must be 
true from (8): 

ail = EUJtJt Y (j) - BUT (9) 


pi) = BUy(®) + BU V ( U > + - 


c T 


( 10 ) 


D = +B 1 ) Ty ( i ; ) +a‘Vi P i ) < n > 

It is interesting to note that the internal dissipation consists of terms 

associated with thermoelastoplastic strains or energy dissipated by plastic 
deformations. 

* . Vc 

To evaluate E'J k * and B'l the von Mises yield criteria and associated 
flow rule may be used as discussed by earlier investigators [10, 11, 121. 
Writing (9) in differential form, 
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where 


daU = Eim(d - d Y <!>) -.BUdT 




( 12 ) 


(13) 


The plastic potential F is related by the equivalant yield stress a and 
the second deviatoric stress invariant J, 

F = a 2 = 3.1 (14) 

from which we derive the incremental yield stress in the form 

■C'-i < 15 > 

and d\ in (13), a positive constanL, can be shown to be 


dX = d v ( v) / 2 a 


(16) 


with the incremental equivalent yield strain dy^ p 3 related by the bilinear 


plastic modulus E as 
<p> 


dv/ = dcr/K 


( p) 


(17) 


Using (13) through (17) in (12) we obtain 

d ^11 = EU^(dY Xj( , - Z kje dY^ p) ) - B^dT (18) 

In view of (15), (17), and (18) d v < v) assumes the form 

dy(p> = If 1 (Z, j^^dy^ - Z^BUdT) (19) 

where 


H- 1 = E + 2 r ,Z ttt Er«*« (20) 

( p ) 

Substituting (20) in (18) yields 

da 13 = EUkjgdy kje + EUkjtdY k£ - BUdT - B»'ldT (21) 


in which 


EUk£ = -ir i E 1 J»»Z Z„ E*.£p<I 

yq B n 

B 1 3 = -ir l B nn Z Z k .E'hi 
mn K Z 


( 22 ) 


( 23 ) 



7 


We now return to (le), a representation of irreversible thermodynamic 
processes and express it for a small time increment in the form, 

p0(At)T](At) - q (At) 3 - ph(At) - D( A t) = 0 (24) 

l J 

Once again dropping (At) and substituting from (10) and (11), it is 
possible to write (24) in the form 

(T n + T)(BU;( + ~) 


qj - ph + - o'ivW = 0 


where 


and 


( P) _ ( « ) 

Yu ~ Yu Y\j 


YU 


Yi j 


y', 5 ’ - BiY 1 ” - Z u ir ‘<Z.,E , "‘Sw - 

Substituting (26) and (27) into (25) yields 

(T 0 + T){(BU +BO + f}U); n + (B + p) T 


cT i * * • 

+ TfT" } - q J - ph 4 EU«» G y. , Y„ 


VC-. , . , x 

E ),n Gv,, v ( B ‘ 

1 1 J 'Bn 


(25) 


(26) 


(27) 


where 


-E Unn W an Yu T+ Ei J ""W mB T Y^* - B 1J Y tJ ‘ P T 


CT 


U(G Y u + Wu T) = 0 


* 


-l * 


pu = h‘buz m z 1 [JI ew^ 


B = H -1 B 1 JZ, ,Z B« 
U pq 

P = -H _1 BUZ n Z pq BM 
* -1 

G = H Z Z E”Pt 

r s pci 

W m = H" l Z Z HPq 

n n m n pq 


(28) 

(29a) 

(29b) 

(29c) 

(29d) 




(29e) 
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Note that all quantities in (28) are now expressed in terms of the total 
strain and elastic strain but not plastic strain. The plastic behavior 

is exhibited by E'3^, B 1J , p 13 , B, p, G, and W nn . It can easily be shown 

i( Vc ~ ~ 

that for isotropic solids B, p, B, B, and are always zero. They need 
be considered only for anisotropic solids. Here elastic strains follow 
simply from the elastic constitutive law. 

The expression in (28) is the transient heat conduction equation with 
a complete elastoplastic coupling and internal dissipation. The formulation 
of incremental finite element equations from (28) will be shown in the 
following section. 

4. HEAT CONDUCTION EQ U AT IONS - F IN IT E ELEMENT FORMULATION 


To introduce the finite element application to (28) the element 
temperature T and element displacements u t (i = 1, 2, 3) are replaced by 
a linear combination of all nodal temperatures T R and all nodal displacements 
u N with suitable interpolation functions [9, 16], in the form 


T = 

(30) 

u i = V uN 

(31) 


For the 8 node isoparametric element, we have a = 8 and n = 24. Q R and 
,|j 1N are the normalized interpolation functions for temperature and displace- 
ments, respectively. 

The strain-displacement relationship is given by 


Yi j = 4 Vi + 

On substituting (31) into (32), we have 


(32) 


A N1] U> 


'NM 1 J 


- YU 


(33) 
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where A and C . are the strain transformation operators, 

N 1 J N M 1 j 

~ ) + ♦jNM) 

^NMl i ~ ^ tkN’ ltkH » J 

Now let it be required to solve the differential equation (28) 
rewritten in the form 

L( YlJ , Yn , T, T) = 0 
or upon substitution of (30) into (33), 

L(uN, un, T R , T R ) = 0 (34) 


where L is the differential operator and L(u N , u N , T R , T R ) is considered 
as the local residual L(u N J u N , i' R > T R ) . Requiring this local residual 
to be orthogonal to the subspace spanned by the functions Q R for each finite 
element j i . e . , 

I 

L(uN, u", T R , T R )n R dv = 0 


I 


(35) 


which is essentially the Galerkin’s method, we obtain the finite element 
model of (28) , 


/ 


{ (T 0 + Q u T^)O r [(BU 4 + pn)(A MU 4- 2 C MP13 u R )u" 


+ CB+B ) n 9 T» 4 “ f3 s T3] “ qj Q„ - php, -I- E ’ J" "Cf} R (A N 


P ( 'e ) 


+ C NMU uN u M)(A Pmrt 4 2 Cpq<A(> u4)uP - E')»G a (A HJ 4 2C pq u t* )u R ^ 


- El^nW B|i n R (A N53 uN + C NMU uV) T t EW..W.AV, 
A n1j uN +2C nmj jU m u n )T 4 0 Tfij, 

" ^ n O R (A NU uN +2C HH1J uNu«)§- = 0 


,( e) 
1 


( 36 ) 


Introducing the linear Fourier law, 

q 1 = 
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where K U is the thermal conductivity, using the Green-Gauss theorem, 
and after some algebra in (36), we finally arrive at the finite element 
heat conduction equations for a kth time increment. 


M T 3 4 . C TS = p(Q) 4 p( Q ) JL . p( C ) 4 p(E p ) 4 p(Tp) 

+ R «» A (k> V) V) V) R (k) 


(k) 

in which 

heat capacity matrix 
conductivity matrix 


N R s j cQ r n s dv 


R S 


volume heat supply vector 


p( Q 

R ( 


/. 

J K lJ n R 


yf 


M n s>J dv 


phO R dv 


(37) 


(38a) 


(38b) 


(38c) 


surface heat flux vector 


>( q ) = / q 1 n 

R( k) J 




(38d) 


Here n 1 is the unit normal to surface. If the surface is uninsulated and 
convection loss is to take place due to ambient temperature T' then the 
following boundary condition should be met: 

qin t + q + cy(T - T') = 0 
or <l 1 n 1 = -q - a(T - T' ) 

where a is the film coefficient. This requires the surface heat flux 


vector (38d) to be replaced by 


P («) 

V) 



O'! ' ) dA 



0 R a?n s dAT s 


(38d- 1) 
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The second term of the left-hand side of (38d-l) may be added to the 
conductivity matrix (38b) so that 


R S 


J H 13 n* M n s> 3 dv + J 


n 9>1 dv + / Of OpfisdA 
A 


and 


p( q ) = - I q 

R (k) J 


(q - o'l")dA 


Pseudo heat capacity vector 


A 


-(If sW - . dv ) T 

( k ) 

Pseudo elastoplastic coupling vector 
p( E p) = ( I (T 0 + T 

R ( k) J (k - l ) 


(k) J A 0 " ° ~ “ 1 ) 

V 


)o„(bu B> ] + I 13 )(a mU 


v 


+2 C u p )dv}u M 

+/ ( k - l) (k - l) 


+ { j( T o H- T ( k _ t) )n R 


(B + fl)n s dv}T3 

(k - 1) 

Pseudo thermoplastic dissipation vector. 


/ 


* * 


P «(k! = ‘ f / E11 '^ (A N>! U ( N k -l 


4- C ,uN U M ) 

X ) N M 1 J ( k - 1 ) ( k - X ) 


(h 0n% + 2C. U<J )dv]u p 

v fnn P5mn (k - i) (k - l) 

t t j b>»CO,(l rll + 20,,,^ _ 


V 


(38b-l) 


(38d-2) 


(38e) 


(38 f ) 


- n 
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f I En^w Q, (A u N -)• C M u |U 'u^)dvlT 
' I b n « N1J NM l J fy.<) tfV 1 (k _ t ^ 

~*v y 

t Jlu.«. A dviT k _ u v< 1 ;> k _ i) 


v 


I 


+2C., .. , ,u M 


) dv] T 


f ( I P n O.R(A Nl ,u" . - 

1 r n (k - x) J (k - l) (k - x) (k - i) 

V 


■/ 


[ l * ^ dV]T ( k- X) + 


/ 


° ,,n » <A ” '>“<"« - n 


V 


+2C U* 

1 J ( k - n ( k - x 


)dv + ( JkV CT 

> Jv“ 




(38g) 


Here (k - X ) refers to the previous time step; and if the pseudo elasto- 
plastic coupling vector and pseudo thermoplastic vector are removed the 
expression (37) is identical to the uncoupled transient heat conduction 
equations. 

5. INCREMENTAL FINITE ELEMENT 
EQUATIONS OF EQUILIBRIUM 


The standard finite element equations of equilibrium is given by 

(39) 


/. 


a 11 i dv = F^» ) 

6u N 


N 


in which F^ a ^ is the nodal applied load vector. The incremental form of 


(39) is obtained by taking a variation of (39) in the form 

f do- 1 J AYu dv + / a i id (AlU. } dv = dF ( a > 

l b u N / ftu N N 

-/V j V 


(AO) 
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The first and second terms in the left-hand side of (AO) correspond, 
respectively, to incremental changes of stresses and geometries. In view 
of (20), (21), (22), and (26) and performing appropriate differentiations, 
it is now possible to write (40) in incremental quantities, 

J. [(E 1 + Eil*A)(A Mltj6 duN + 2C NMkje uMdu H ) 

- (BU + b 1J )n R dT R ] (A N1 , 4 2C MHU u")dv 


l 


2 ct U c mi1 , du M dv = dF (a) 


NM j j 


(41) 


It should be noted that CT iJ in the last term of (41) implies the initial 
stresses or the residual stresses in the structure just prior to a new change 
in geometry. 

After some algebra the final form of incremental thermoelastoplastic 

equation of equilibrium for the jth incremental step 

(K ( « ^ + k ( G ^ 4 K ( p ^ )du M = dF ( a ^ + dF ( T ) + dF ( n) 

. K NM ^ ^NM ^ N (J) N ( J) N (J) ' ^ 

in which K^, K ( N «>, and K^p> are the standard stiffness matrices represent- 
ing linear elastic, geometrically nonlinear, and plastic behavior, respectively. 


K { « 1 = 

NM 


f mill* 

J v 

(g) = r 

NM I 

J v 

-/*■ 


2 ° UC N M 1 J dV 


J kf'A A , /i dv 

n Ni j Mkf v 

and the incremental thermoelastoplastic load vector dF^ 1 ^ is 

<;’=( / C" + » U )n»(A„u * 2C„ n u» )dvldT« 

( j ) / ( i - i ) ( j - % ) 

J V 

+ (f B1J 0,C NHU dv}duM T* 

I ( j - l ) ( j - X ) 

•2 v 

All the rest of the terms in (41) other than those mentioned above may be 
grouped in F (N ^ called the pseudo nonlinear, load vector but may be dropped 

V ) 


(A3) 
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because of their negligible effects. 


6. SOLUTION PROCEDURE 

We are now provided with incremental heat conduction equations and 
equilibrium equations. Either thermal loads or mechanical loads or 
both may be applied. Depending on loading and boundary conditions, we 
can either start from equilibrium equations or heat conduction equations, 
but both equations must be solved iteratively within a time increment. 

Any existing recursive formula for step-by-step integration or difference 
operator may be used to solve heat conduction equations. In the present 
study a difference operator for linear variation of temperature within 
a time increment is combined with iterative solution of plastic equilibrium 
equations. 

The incremental temperature at any time step k is given by [14,15]. 


T = 2 (~ N + 8)-i [P 

At -- ~ - ( k 
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+ -ST } - T (44) 

At — ( k - X ) k - l) 

, / E p ^ 

Here N, £, P, and T are assembled forms of N nm , K nm , Pj^ + P ( R q ^ + P r + p r 
+ PjJ rp ) and T R , respectively. The reference temperature and initial thermal 
input can easily be incorporated in (44) and the temperature change at the 
end of the first time increment calculated. 

The results of this solution are used in the assembled incremental 
equations of equilibrium to determine the displacements and stresses. 

These stresses are checked, element by element, for yielding. If any 
element has yielded the plastic tangent stiffness matrix is constructed 
and the standard iterative cycles are repeated until convergence is achieved 
[10 , 11 , 12 ]. 
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With the final values of displacements, it is possible to calculate 
the displacement rates by 

u it = (u k - u k -3.) /At (45) 

for use in the heat conduction equations, Returning to the heat conduction 
equations for the second time increment, the process is repeated as before 
except that the elastoplastic coupling is now to be included. The elasto- 
plastic and thermoplastic model as determined in the converged solution of 
the plastic equilibrium equations of (42) will be used in the heat con- 
duction equations. The marching with time increments, thus, continues 
until the desired length of time has been reached. Details of the 
computer program are given in Appendix A - Capabilities and Limitations 
of the Program, Appendix C-Flow Chart, Appendix D - Subroutines Organiza- 
tion Chart, Appendix E - Description of Subroutines, Appendix F- Data 
Input Format, and Appendix G - Program Listing. 

In the present study, we use temperature and displacement approxima- 
tions based on a three-dimensional linear isoparametric function and the 
integration is performed by an 8 point Gaussian quadrature [16] (see 
Appendix B) . 

7. APPLICATIONS 

In order to verify, first of all, correctness of the present 
approach a comparison study was made for an uncoupled heat conduction 
of a beam reported by Wah [17] who used a classical series solution and 
substantiated his results with Boley [19]. An excellent agreement was 
obtained as shown in Figure 1. 
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Next, a series of example problems were tested to determine effects 
of various terms in the governing equations; namely, behavior due to 
coupling and non-coupling, linear elastic and elastoplastic properties. 
Figure 2 shows the information on geometry, boundary conditions, material 
constants and temperature input. Only one-quarter of the symmetrical 
three-dimensional solid is shown. The temperature change of 200°C is 
applied at nodes of the center-left end element. The material properties 
given in Figure 2 represent a mild steel. The transient temperature dis- 
tribution in the direction of x with y = z = 0 is shown in Figure 3. 
Effects of elastic and elastoplastic couplings are studied in this case. 

It is interesting to note that there exists little difference in tempera- 
ture distribution between elastic of elastoplastic coupling for the time 
period examined. However, the displacement w at y = 100 mm and z = 300 mm 
becomes larger for elastoplastic coupling than that for elastic coupling 
after approximately one hour as shown in Figure 4. The changes of w dis- 
placements vs. time at point A are plotted in Figure 5. It is seen that 
uncoupled displacements are larger than the coupled displacements, a 
trend confirmed by Oden and Poe [18] in their study of thermoelastic 
one— dimensional problems. It should be noted that the elastoplastic 
displacements are larger than the elastic displacements in both cases. 
Figure 6 reveals an interesting fact that temperature is also lower for 
a coupled case than for an uncoupled case but little difference is noted 
for either elastic or elastoplastic behavior. Figure 7 shows the stress 
in z direction o z from that for elastic coupling, a fact well known in 
mechanics. Finally, plastic regions developing with elapse of time are 
shown in Figure 8, indicating little effects of coupling. 
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The effects of convection through uninsulated surfaces with the film 
coefficient a - 1.0 Kg/mm hr °C, ambient temperature T' = 1000°C and the 
heat flux q = -100 Kg/mm hr at z = 300 mm on the x-y plane are investi- 
gated and the results are shown in Figure 9. Large elastoplastic 
displacements (w) result due to surface heating together with ambient 
conditions. Variations of material properties from element to element 
are accommodated in the program and an example problem for such case 
is described in Figure 10. Once again the transient temperature dis- 
tributions are almost identical for elastic coupling and elastoplastic 
coupling as shown in Figure 11. Significant deviations exist, however, 
for displacements (w) (Figure 12) between elastic aud elastoplastic 
couplings as temperature rises as noted earlier for , the uniform material . 
Because of ambient temperature and heat flux through uninsulated surfaces 
and variable material properties throughout the structure, the pattern 
of development of plastic regions (Figure 13) differs considerably from 
that of the previous example of uniform material and insulated surface. 

In the foregoing example problems, geometric nonlinearities are excluded 
for the interest of computing time. 

It should be noted that for the case of an isotropic solid, the 
tensor of thermoplastic moduli (23) and expressions of (29a), (29b), 

(29c), and (29e) are zero lint must be updated in the case of an aniso- 
tropic solid as plastic deformation progresses [9]. Although the 
temperature dependent material properties, finite strains, and the 
dynamic-coupled inertia effects can easily be handled, such applications 
are not included in the present study. 
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DISPLACEMENT BOUNDARY CONDITIONS : 
u «v-w=*0 At x»0 and x-600 mm 
v-0 At y-0 

w-0 At z-0 

TEMPERATURE BOUNDARY CONDITIONS: 

Insulated on all the surfaces except at 
the points (x,y,z) = (0,0,0), (0,100,0), 

(0,0,100) and (0,100,100) which are 
kept at 200°c. 

CONSTANTS ; 

E = 2.0 x 10* (kg/ram® ) , E = 2.0 x 10 3 (kg/mrn* ) , o Y * 25 . 0(kg/mm») 

(p) 

v “ 0.3 , a = 1.3 x 10 _f3 (/°C) , K - 9.0 x 10 3 (kg/hr °C) 

c - 0.3 (kg /mm* °c) , T 0 = 27(°C) , At =0.05 (hrs.) 

Figure 2 Discretized geometry of three dimensional solid and input data. 
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Figure 9: Comparison of Coupled Elastoplastic Displacements (w) With and Without Surface 

Insulation at y = 100 mm, z = 300 mm in the k-Direction. 
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DISPLACEMENT BOUNDARY CONDITIONS : 

u = v = w = 0 At x = 0 and x = 600 mm 

v = 0 At y = 0 

w = 0 At z = 0 


TEMPERATURE BOUNDARY CONDITIONS : 

T = 200°C at the points (x,y,z) * 
and (0,100,100). 

a = 1 . 0 (kg/mm. hr °C) , T* = 1000.0 
on the sufface 7. = 300 mm. 
Insulated on all other surfaces. 

CONSTANTS : 

E = 2.0 x 10 4 (kg/mm 2 ) 

E = 1.0 x 10 4 (kg/mm 3 ) 

E = 0.7 x 10 4 (kg/mm 3 ) 

E 


( 0 , 0 , 0 ), ( 0 , 100 , 0 ), (0 0 , 100 ) 

(°C) and q = -100.0(kg/mm.hr) 

for elements 1 ~ 6, 
for elements 7 ~ 12, 
for elements 13 18, 

= 9.5(kg/mm 2 ) , v * 0.3 

(kg/hr °C) 

At « 0.05 (hrs . ) 


( ?) 


= 1.0 x 10 3 (kg/mm?) 


o\ 


a = 1.3 x 10-s (/ C) 
c 


9.0 x 10 3 


0.3 (kg/mm 3 °C) , T 0 = 27 (°C) - , 


Figure 10: Discretized Geometry of Three-Dimensional Solid and 

Input Data. 
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Figure 11: Temperature Distribution of y = z = 0 in Figure 10. 
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8. CONCLUDING REMARKS 

A three-dimensional thermoelastoplastic analysis has been carried 
out using the incremental theory consistent with the first and second 
laws of thermodynamics. Complicated mathematical operations emanating 
from the functional theory or state variables are replaced by proposing 
an incremental free energy as a function of total and inelastic strain 
and temperature unique only within a small time increment. The similar 
incremental functional dependency is then valid for stresses, entropy, 
and heat flux. Such treatment lends itself to numerical techniques 
taking advantage of the finite element method and time integration by 
difference operators. 

For the example problems and material properties considered in this 
study it appears that elastoplastic coupling is significant for dis- 
placements and stresses, but that neither elastic nor elastoplastic 
coupling has any effect on temperature distribution. It may be argued, 
however, that for certain types of material and geometry these conclusions 
would not necessarily be true. Exhaustive study on such effects is 
beyond the scope of the present report. 
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APPENDIX A 

CAPABILITIES AND LIMITATION OF THE PROGRAM 

General. This program analyzes a three-dimensional solid subjected to 
both thermal and mechanical loadings. The program takes into account 
the elastoplastic behavior coupled with transient heat conduction. The 
formulation is based on the first and second laws of thermodynamics, von 
Mises yield criteria, associated Prandtle-Reuss flow rule, and the linear 
Fourier law. The finite element discretization by means of a linear iso- 
parametric interpolation for both temperature and displacement fields is 
utilized. Integration is performed by Gaussian quadrature. A step-by-step 
time integration assuming a linear variation of temperature within a time 
increment is used to solve heat conduction equations. Capabilities and 
limitations of this program are listed as follows: 

(1) Capable of handling up to 60 nodes and 30 elements. 

(2) Temperatures may be specified at nodes (50) . 

(3) Heat flux and heat supply may be specified on the element surface 
(50) and inside the entire element solid (30), respectively. 

(4) Surfaces may be insulated or exposed to ambient temperatures. 

(5) Capable of incorporating 100 restrained generalized coordinates. 

(6) Geometric nonlinearities are not considered in the program. 

(7) Capable of handling laminated structure with varying material 
properties from element to element. 
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APPENDIX B 

VARIOUS INTEGRALS IN ISOPARAMETRIC ELEMENT 

In the present study, linear hexahedral isoparametric elements are 
used to model the three dimensional solids and consequently constitute 
the basis for displacement and temperature fields. Although details of 
isoparametric elements may be found in Zienkiewicz (1972), some of the 
integrals essential in the program by Gaussian quadrature are shown in 
explicit form: 

/v'|fN d ' / = i a f* (1 4 SSnHI + Tm M )(l + ??M)|j|d?dT)dC 

/t* *tn dv = f t f f % (1 + ??„)(! + Tm N )d + CChH 1 + + TTHh) 

( J + CC M )| J l d ?, dT 1 d C 


/ > y'|/ N 'llM'l(R dv = "g 3 /; /; r (1 


+ SSmH 1 + N ) ( 1 + CCn) (1 + ^ M ) 


(1 + T|T)m ) ( 1 + C£m)(1 4 + 31T] R ) ( 1 + QQ n ) 


X J| d^dTjdf 


./v* NAwijdv = | // /; /; O + SSnHI + Tin N > ( 1 + CCN)Al,.i J | j|d?dT)d C 
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Jv' |fN'l'q A MU dv = I s i: i: i: (1 + SSnH 1 + T)T1 N ) C 1 + CCn) 

(1 + SS Q )(1 + T]T1 q )(1 + CCq) X A Mlj |j| d^dTldc 

/»♦*, l*H,J dv '“ // // J/ , J i J l d ? d< n d c 

/v A HBn A NU dv = £ £ A Mon A NU |j| d ? dTld C 

/a N dA = | /; r C 1 + ^nH 1 + T)T! N ) ( 1 + £Cn) dA 

and 

/a iMM dA “I® ^ / 1 ^ + ! ^ n ^ 1 + + CCn) 

(1 + ?? H )(1 + HT1 m )(1 + CCm) dA 


where the differential area dA is expressed as 



g 9 

MM , (M M M , fox by W 

bT) b? biy + \b? bTi b? bTi; \b| ?>ti be biij 


k 

de;dT| 


for Q = ± 1 plane 


r bZ M M .M 3 , (hi M M 2 , /bx M M MV 3 

>.m) sc " ^ bcj v*>*n be " sir &d \»n »c “ .wi &c/ 


* 


dT|d C 


for 5 =i 1 plane 
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/JiY „ & £l\ 

2 Jte M 

. £x zA 

8 

.■fix £xf 

\tQ 6? 65/ 

V>C &S 

bc ?>?/ 

V&C B? 

Be B?y 


for T] *» ± 1 plane 


Here, displacement and temperature fields are related by 

8 8 
u = £ vj; ! u , etc., and T = E i|f i T t 
i=l i— 1 


h 

d C d S 


where 

- I (1 + SSi) (1 + Tl^!) (1 + CCi > 


The determinant of Jacobian is given by 
8 8 8 

l J 1=7P s, .S, ,S, pM 1 + TITliHl + CCiK CTlj(l + ??j)(1 + CCj) 

8 1=1 j=l k=l L 

Yj&c(l + ?? k )(! + TlTl^Z, 

- tm 1 + sSjM 1 + cc^ZjcX 1 +^ k >( 1 +^k)yk] 

- ? t (i + T|T) i ) ( i + cCiXitV 1 + ss 3 > (1 + cc^XjCkC 1 + ss^x 1 + 

- 71,(1 + + CCpZjCkd + (SkX 1 + WicK} 

+ ? t (i + jiti t > ( i + cci)MVi + «,xi + ccj ) x jCk (1 + ??,)(! irn k )yic 

- TlX 1 + SSjH 1 + CCjXjCX 1 + WkX 1 + T 1V x k}] 
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APPENDIX C: FLOW CHART 


INITIALIZATION 


Time step (k) = (), Temp^ = <>, <5^^ - Yield Strees(a Y ), u ( k =0) 


k -t ! | 

Calculate T^ from Eq.(35) 


Calculate u ^ from Eq. (32) 
Calculate strains and stresses 
Calculate equivalent stresses, a. 


; ~ X 

^(k) ■"* CT max X 


( 1 ) _ 

da (k) ®(k) ~ a (k-i) 


Calculate ft ^ from Eq. (27) 
and pseudo plastic load dF^ from Eq. (33c) 


Calculate elastoplasti'c coupling vector @ k - 2 

i 


Calculate from Eq. (33) 


dT (k) “ - f (k) ' ~(k-l) 


Calculate thermal load vector dF^^(Eq. (34)) 


du., , = K _1 fdF^ P ^ + dF (T) } 

-(k) - - ( k ) 
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APPENDIX D 


SUBROUTINE ORGANIZATION CHART 


MAIN" 


DATA 


-Cl 


ZERO 


MCODE 


PRINTT 


P ST IFF 


ASSIGN 


REDUCE 


DFVCTR 


-[matmpy 


DTHETA 


RSTORl 


PRNTDT 


STRSEL 


- skappa] - 


h 1 STIFF | 


ELASTC 


7ER0 




ZERO 


— |_ CONTH ; 

— |;;mcodf; 


ZERO 


K STACK 


-| PR INTK | 


MATINV 


_ L ASSIGN | 


-f MA TMPY 


REDUCE 


ZERO 


RESTOR lj 


— D^ ro J 

— f ZERO 1 


r 


- ['GAUSS | f F 


PLAMDA 


PMU 


PHI 


- PJCOB 


i — 1 PLAMD A 1 


'HST RES 

ZERO 






"1 PMU 


nATKS 

If | 


r : r 



i 1 


PHI 


PJCOB 


STRSPL 


ZERO 


PHI 


PMU 


PLAMDA 


PJCOB 






















41 


Subroutine 

Name 

ASSIGN 

ROUND 

CONTH 

DATA 

DFVCTR 


DTHETA 

ELASTC 

F 

GAUSS 

HSTRESS 

MATINV 

MATMPY 


APPENDIX K 

DESCRIPTIONS OF SUBROUTINES 
Descriptions 

Rearranges the nodal displacement vector 

Applies the boundary conditions and reorders 
the matrices accordingly 

Defines various quantities for Guassian 
Quadrature integrations 

Reads all input data 

Calculates all psuedo coupling vectors and 
heat input vectors 

Solves heat conduction equations 

Calculates the elasticity matrix 

Function subroutine for Gaussian integration 

Integration by Gaussian quadrature 

Calculates the thermal load vector in 
equilibrium equations 

Matrix inversion 

Matrix multiplications 
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Subroutine 

Name 


MCODE 


K STACK 

PHI 

PJCOB 

PLAMDA 

PMU 

PQR 

PRINTK 

PRINTT 

PRNTDT 

PST IFF 
REDUCE 

RESTOR 1 


Descriptions 

Reassigns the global node number to the local 
element node number 

Assembles the local element matrices into a 
global form 

Function subroutine involved in Gaussian 
integrations 

Calculates Jacobian in Gaussian integration 

Function subroutine involved in Gaussian integration 

Function subroutine involved in Gaussian integration 

Defines constants necessary for Gaussian integration 

Print stiffness matrix, heat capacity matrix, and 
conductivity matrix 

Print nodal temperatures 

Print nodal displacement and equivalent nodal 
forces 

Calculates the plastic tangent stiffness matrix 

Modifies the equivalent nodal vectors in 
correspondence to BOUND 

Restores the equivalent nodal vectors to the 


form prior to REDUCE 
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Subroutine 

Name 

SKAPPA 

STIFF 

STRSEL 

STRSPL 


Descriptions 


Calculates 


Calculates 


Calculates 


Calculates 


the plasticity matrix 
the elastic stiffness matrix 
stresses prior to yielding 
incremental stresses after yielding 


ZERO 


Zeroes out all matrices 
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Card 


Card 


Card 


APPENDIX F 
DATA INPUT FORMAT 


1: 

FORMAT (7 15) 



(1) 

NELEMT 

- 

Number of elements 

(2) 

INODE 

- 

Number of nodes 

(3) 

NB 

- 

Number of constrained displacements 

(4) 

IPT 

- 

Number of integration points in the 




Gaussian quadrature 

(5) 

NBHC 

- 

Number of prescribed nodal temperatures 

(6) 

ITER 

- 

Number of time increments 

(7) 

NKE 

- 

Number of type of elements 

2: 

FORMAT (3 F 10.0) 



(1) 

RT 

- 

Reference temperature 

(2) 

EPSS 

- 

Percent error limit, e in the elastopla; 




analysis 

(3) 

DELT 

- 

Incremental time interval, At 

3: 

FORMAT (8 F 10.0) 



(1) 

!■:(!) 

- 

Young's modulus 

(2) 

EP(I) 

- 

Plastic modulus 

(3) 

SYIELD (I) 

- 

Yield stress 

(A) 

XNU(I) 

- 

Poisson's ratio 

(5) 

TKX(I) 

- 

Thermal conductivity 

(6) 

SH(I) 

- 

Specific heat 

(7) 

ALPHA (I) 

- 

Coefficient of thermal expansion 

(8) 

DENSTY (I) 

- 

Density 



45 


Cards 4: FORMAT (16 15) 

(1) NY (I) - Type of elements I 

Repeat cards 4 as required to complete all elements. 


Cards 5: FORMAT (815) 

(1) to (8) 

MA(I), MB ( I ) , MS ( I ) - Node numbers of element I 
Repeat cards 5 NELEMT times. 

Cards 6: FORMAT (6F10.0) 

(1) to (3) 

X(I), Y (I) , Z (I) - x, y and z coordinates of node number I. 
Repeat cards 6 as required to complete all nodes. 


Card 7: FORMAT (2 15) 

(1) LSTRES = 0 

= 1 

(2) LHEAT = 0 

= 1 


if stress anslysis is not desired (no elastic 
or elastoplastic coupling) 

if stress analysis is to be included (elastic 
or elastoplastic coupling present) 
if heat conduction analysis is not desired 
otherwise 


Card 8: FORMAT (2 15) 

(1) KLOAD = 0 

= 1 

(2) MLOAD = 0 


if mechanical loads are not applied 
if mechanical loads are applied 
if temperatures are not prescribed 
if temperatures are prescribed 


1 



(An assumption is made LSTRES = LHEAT = MLOAD = 1 and KLOAD - 0 
after card 9) 

Cards 9: FORMAT (15, FlO.O) 

(1) ILIBND(I) - Node number at which temperature is prescribed 

(2) TFR(I) . - Prescribed temperature 

Repeat cards 9 NR1IC times. 

Cards 10: FORMAT (4l5, 3F10.0) 

(1) * NS(L,I) - Node numbers of the uninsulated surface (I 

changes from 1 to 4) 

(2) * SC (L, I) SC (L, 1) = o' (film coefficient) 

SC(L,2) = q (heat flux) 

SC(L,3) = T 1 (ambient temperature) 

*L indicates the number of uninsulated surfaces subjected to convection. 
Provide a blank card to signify end of data. 

Cards 11: FORMAT (4l5) 

(1) NODE - Node number at which displacements are constrained. 

(2) IU =0 if U is not constrained 

=1 if U is constrained 

(3) IV =0 if V is not constrained 

=1 if V is constrained 

(4) IW =0 if W is not constrained 

=1 if W is constrained 

Provide a blank card to signify end of data. 



47 


ards 12: FORMAT (8F10.0) 

(1) D11V(I) - Heat supply of element I. 

Repeat cards 12 as required to complete all elements. 



APPENDIX G 


COMPUTER PROGRAM LISTING 


COMMON JbLKA/MELEmT, INODE, Nb.NBHC,LHEAT,LSTRES,nN,NM,NNH,NMH 

COMMON ;bL)CR/N0PRNT,NPa6E,LINES,MN,K0DE.KT0TAL,mT0TAL,NL0AD,KL0AD 



TYPE Jt 8 ,8 ) ,TYPEM I 8 , 8 ) «TTPEN I 8 ,8 ) .TYPES l B ,8 ) ... 

COHMON/bLO/SK 1200001. HE 12000 1, CKI 20001. HCKI 2000). 1BNDI 100). 
|HBNO<SO> 


7 * 


B» 

»• 

10 * 


IIP 

I3» 

IB* 

IS* 

u* 

I7» 


IB* 

20* 

21 * 

22 * 

2 J» 

2 B» 

25* 

»*• 

27* 

2B» 

n* 

» 0 » 

21 * 

»• 

»• 

H* 

»• 

)l» 

37» 

36 * 

»• 

BO* 


B|P 

B2» 

B»p 

BB» 

B 5 » 

B*» 


<I7p 

BB* 

BB» 

SO* 

S|P 

52* 

SJ» 

S»P 

$ 5 * 

SB* 


57* 


58 * 

SB* 

40 » 


*!• 

* 2 * 

* 3 * 


COHMON/rLK/NAI 30) ,HB I 30 > , MC I 30 > ,H0 1 30) ,MP I 30 ) » Mj(S0l t MR (30 1 
COMMON /BLKI/AI 101,81 10), C(10I,PIS3I,9|S3),R(3S).0|*).MI*I, 
COMMON /BLK2/MII8 | ,M2I8) .HAIR) ,MR(B| ,HSI8I ,H4I8) ,H7|8) ,H*I8 
pm? 1 8 1 ,H|0( 8 | ,61 (8» ,S2(B) ,63181 ,6B|8l ,65(8) ,6818) ,S7IB) ,«BtB 
COMHON^bUKM/TOI « 0) , Til 40) ,721801 ,DT0t40> ,0T1 (AO) ,0721401 
C0HM0N/;LK3/NS|S0,SI , SC 1 80, 3) « AREA BIBO) 

C0HH0N/8LICS/t6AUS,ITER,EPSS»1SE ...... 

COMM0N/bLK 6/0(4,8) ,Et 10) ,EPI 10) .STIELDI IO|,XMUI|0),TE2IIO), 

• SHI 10) .ALPHA I 10) .DEMSTYI 10) , NY I 30) ,DELT ,MT 
C0MM0N/bLK7/X(60) ,Y(40) ,2(40) .MNODEIB) ,MK0DEI2B| 
C0MMON/bLK8/00( 1 80 ) ,Ul 1 1 80) ,U2 I 1 801 , OUO *• 80 1 YOU I II 80 I .01121 I 
COMMON /BLKV/DS1SHAI4) ,SIGHAI4,30> iSISMAl 14,30) 

C0MM0N/BLK i6/SISMAB(30) .SISMIBISO) ,DS1 S 1 B 1 30) ,Ds I 62B I 301 , 

• S 1 GN A X I 30 ) . 

COMMON/bLK) |/STRN|A, 30) »DSTRNIA,30) ,5TRNPI4,30) ,ABC 
DIMENSION Cl I I 80 » ,BSK I 1 80) ,07140) , DEB (40) ,0HVI30> .EOF I 8) 

01 MENS | gN OP I 2B ) , FR I I 80 ) , PR • 2B , 2B ) , A Jl 4 , 4 ) ,FP 1 1 80 ) , TP I 4Q ) 
lYIELD.p 

TIME-0. 0 
II 1-0 
NNN-0 

CALL DAT*IFP.BSK,C1 I 
00 IS i.i.nelemt 
NY 1-NYI j> 

S1SMAX|j)mSY1ELDINY1) 

is continue 
ne-nELEmT 

DO 12 |.I,NNH 
TPI l )PTpl I ) 

T2U)*Tg(ll 
12’ Til I )PTOl I ) 

CpppppTIHE INCREMENT 
1000 UI-MIpI 
2000 NNN»NNn»I 

TIMEpTjmEpoELT 
IFIlII.gT.lTERI STOP 
■R 1 TE! A, A ) 111, TIME 
REAOIS,*) (DHVMI.I-I.NELEHT) 

13 CONTINUE 

RRITEI4.10) (0HVII1,1p|,NELEMT) 

DUN-0. 

CALL 2ER0I0F', INOOE.I) 

L-l 

NX-0 

20 CONTINUE 

CALL MCpOElMKOOE,MNODE,L,3l 
DHpOHV(l) 

ICAPpI, 

CALL OFyeTR(EOr,RAP,OM,NX,NNN,l,Ll 
00 30 1-1 ,8 
HNO-MNOpEI I ) 

30 0FIMN01.E0FI 1 IpOF(MNO) 

IFIL.GE.NElEHT) 60 TO <10 

L*Lp 1 

Nx-i 

IFINTlL)*NE.NY(L-i)l NX»0 
60 TO 20 


.METSO) 
IPT 
) , 

) 


801 





69* 


45* 
66 • 
6 7* 
6 B • 
61 • 
70* 
71* 
72* 
73* 
71* 
75* 
76* 
77* 
78* 
79* 
80* 
81* 
82* 
83* 
8 1 * 
85* 
86 * 

• 7* 

• 8* 
89* 
90* 
91* 
92* 
93* 
91* 
95* 
96* 
97* 
98* 
99* 



90 CONTINUE 

c»***«surmce heat boundary 

L * I 

91 IF (NS (L I I) ,EQ*0> 60 TO 92 
AI*AREA9<L)*(SCfl,2)-SC<L,l )*5C<L,3) I 
DO 93 |-1 ,M 

NT»NS(l", 1 » 

93 OMNI 1*oF (N! 1-A1 
L *L ♦ l 
GO TO 91 

92 CONTINUE 

WRITEU,*) , I ,DF( D • 1*1 • INOOt) 

C*****OTHET A yEC TOR 

CALL OThETaIOFB.OF ,MLOaO,TO»TI ,TP> 

CALL PRtNTt(DF,DFB.X,Y,Z,1 , I NODE , NP AGE , 60 ) 

C*****OT VECTOR 
NX-0 

CALL ZEROIFR. 180. I > 

IFINY ( li .N e»NY ( I - 1 ) 1 NX*0 
NYl-NY( j ) 

ALPHA I -ALPHA (NY I >*E(NYI J/C 1 ,-2.*XNUfNY| ) > 

CALL MCoOE 1 MKOOE .HNODE , I ,3) 

CALL MSTRES C TO. FR.MNOOE.MKODE, NX, ALPHA I ,60,180) 

NX* 1 

260 CONTINUE 

DO 265 i*| « NN 

265 FR1I )-FR( I »*FP< ! )*FL0ATI 1 t | ) 

1 F ( IGE.eO'O) go To 269 

C*****GE0METr JCALLT NONLINEAR VECTOR 
CALL ZEpO(AJ.6,6> 

DO 264 j-I.NELEMT 

CALL HC00E<NK0DE.MN00E,I ,3) 

CALL ZER0(Df»,29,l I 

CALL PSTlFF(PK,AJ,t ,IGEl 

CALL AS5!GN<PK,NN0DE.0,8,! t 29> 

DO 267 j- I , 29 
00 267 i( ■ I ,29 

NF-MKODE <F) 

DUO I *U | f NK ) 

267 0P(J)*0P( J)*PK( J.K)*PU01 
DO 268 j* I ,29 
NK«MK0DE<J> 

260 FRlNK)-oP( J)*FR(N k ) 

266 CONTINUE 
269 CONTINUE 

C • • • • *DU VECTOR 

CALL REdUCe(PR. IBN0,NN,N8,|80) 

CALL NatRPT(SK,NN,FR. U0,I80> 

CALL RSTORI IFR, IBN0,NN,NB, 180) 

CALL R5tOR|( UO. I0N0,NN»NB, 180) 

CALL PRNTOTCFR, Uo , X , Y , Z , MM , 1 NODE , NP AGE . I 80 ) 

C*****STR A I N /, STRESS 
WRITE (6^88 ) 

DO SO i.i . neLemt 

CALL McoOC(MKODE#MNODE,I *3) 

CALL STRSELl UO. TO. I YIELD, I ) 

RR I TE ( 6 , 87 ) I , I Y lELD.SlGMAB ( I ) , (SIGMA ( J, 1 I , J*l ,A» ,SI6HAX|I» 


121* 
122* 
123* 
129* 
125* 
|26* 
127* 
128* 
J29* 
130* 
131* 
I 32* 
133* 
I 39* 
1 35* 
I 36* 
137* 
J38* 
1 39* 
1 90* 
| 9 I • 
I 92* 
193* 
199* 
1 95* 
(96* 
1 9 7* 
I 98* 
199* 
150* 
1SI* 
152* 
153* 
159* 
i55* 
156* 
IS7* 
158* 

I 69 * 
160* 
i 6 1 • 
162* 
|63* 
J69* 
165* 
166* 
167* 
|68* 
J69* 
170* 
171* 
172* 
A 73* 
1 79* 
1 7S* 
176* 
177* 




50 CONTINUE 

IF( I YIElD.EQ. 1 ) GO TO 3000 
99 CONTINUE 
C • • • • • PR I NT 0|>T 
WR I TE ( 6 ^ 7 ) 

DO 120 1*1, INODE 
J-3* 1-2 
JI*J* I 
J2*J*2 

120 WRITE (6^66) I ,U0< J) .U0( Jl ) .U0( J2 ) ,TO(I ) 

C • •• • *SH I F T 

00 125 j*|, INODE 
T 2 ( I)*T|(I) 

T 1 ( I ) *Tp ( I ) 

DO 125 j*l,3 
I J* ( I • 1 ) *3* J 
U 2 ( I J ) ■ u 1 ( I J ) 

U1 ( I J) -uO( | J) 

125 CONTINUE 

DO 126 J * 1 , NELEMT 
S1GMIBI i )-SlGMAB( I ) 

00 126 j-|,6 

SIGMAI | j , I J.SIGMA1 J,I ) 

126 CONTINUE 
GO TO loOO 

C • ••• *F OR V I ElDEd 
3000 CONTINUE 
NNN-0 
0SIG2-0, 

00 130 !*|, NELEMT 
0SIG2BC i )-•! ,0 

IF (SI GMaB( I) • GE • S I GM A X ( I ) ) DS I G2B ( I)*I ,0 
DSIGIBf f ) *S I GM AB ( I ) -5 1 GM I B ( I ) 

130 DSIG2*DS1G2*DS!GIB( !)**2 
3500 CONTINUE 

nnn-nnn# I 

IF(NNN-gT.To) stop 

WR I TE ( 6^ 9 ) I I I , NNN 
CALL SkaPPA 
REWIND 2 
KAP-I.o 
9000 CONTINUE 

CALL ZERO I OF ,60 . 1 ) 

REWIN02 
L« I 
N X *0 

25 CONTINUE 

CALL MCqDE(HK00E.MN0DE,L.3) 
oh*dhvi l i 

CALL DFvCTRCEOF.IfAP.DH.NX.NNN.O.LI 
00 35 1.1,8 
mno-hnooM I ) 

35 DF(MNO).EDF( I )*DF(MNO) 

IF(L.GE.NClEMT) Go TO-95 
L*L* I 

IF(NY(li .NE.NY(L-I ) ) N X * 0 



178* 


60 TO 25 


|79» 

Ieo» 

i e i • 
i«2* 

1»3» 

|87» 

1*5* 

?84* 

I«7» 

166* 

!*»• 

i*0» 

1*1* 

I 72* 
193* 
175* 
J»S». 
174® 
J| 7 7 • 
178* 

I 77* 
200 * 
201* 
202 * 
203» 
207» 
205* 
20*« 
2077 
20*7 

2077 

210* 

211* 

2II» 

2IJ7 

2177 

*15* 

21 ** 

217* 

2I»7 

217* 

2 * 0 * 

2*l» 

***• 

2 * 1 * 

2 *«° 

2*5* 

2*5° 

2*7° 

**•• 

2*7» 

230® 

**1* 


2*7® 


75 CONTINUE 

WRITE<6,SI I I ,DFl I ) , !•! • tNODEI 
{•••••SURFACE HEAT BOUNOARY 
L*l 

74 IFINSlLi I • .EO.OI GO TO 77 

A I -AREA. 7It»»ISC I L, 2) -SCIU, I I *SC IL , 3) I 

DO 78 1.1 ,7 

N I “NS (L j I ) 

7B OMNI I.oFINI l-AI 
L«L*I 
GO TO 74 
77 CONTINUE 
{•••••OTHET A yECTOR 

CALL OTmETaIDFB.OF.MLOAO.To.TI ,TP) 

CALL PR i NT T I OF , DFB , X , Y , Z , | , 1 NODE ,NP AGE . 40 ) 

VECTOR 

00 250 j»l, INODE 
OTO I I l«TO( 1 I-TI II ) 

250 CONTINUE 

CALL ZErOIFR. 180. I ) 

00 241 i*l, NELEMT 
IF(NY(Il.NE,NY(l-I»l N*«0 
NVl*NYlil 

ALPHA|.»LPhAINYI UEINYI I/I I «-2,*XNU|NYI I ) 

CALL McjjOElNKODE,NNODE,I,3l 

CALL HSTRES IOTOiFR ."NODE iHKOOE .NX , ALPHA I ,40, I Bnl 
241 CONTINUE 

00 270 I • I , NN 
270 FRI I >*FR< I «*FP< I I 
(•••••OP VECTOR 
RERIND j 
RE*INO 5 

DO 1 74 t*l,NELE"T 
RE AO i 21 AJ 

CALL HCOOE I HKODC , HNODE , I , 3 I 

CALL ZERO! OP ,27,1 I 

CALL PSTIFFIPK.AJ, I, IGE1 

CALL AS;tGN(PK,MN00E,0.S,t,27l 

00 175 J*I ,27 

00 175 K«l ,27 

nk*hkode<ki 

DUO I "UO I NK I “U 1 INK | 

IFINNN.f 8.II OUOMUI (NKI*U*INK» 

175 0P( JI«DPC J)»PKI J.KHOUOl 
00 177 J* 1,27 
NK*MKOOE ! J I 

177 FRINKI.oPI J)*FR!NKl 
I 74 CONTINUE 
(•••••DU VECTpR 

CALL RE0UCEtFR.«8N0,NN,NB, I SO > 

CALL HAT«PYISK,NH,FR,OUO, IB01 
CALL RStORIIFR.IBNO.NN.NB.IBO) 

CALL RsfORilOUO,tBNO,NN,NB,I*OI 
CAUL PRhTOTIFR,DUo,X,T,Z,NH, I NODE, NP AGE. I 80 1 
(•o*«*STRA IN 5 STRESS 
RERIHO 2 


235* DO 170 j*l, NELEMT 

234* CALL HCODE(HKOOE,MNODE,I,3I 

237* RE AD I 2 1 AJ 

238* CALL STrSPl(DUCi,DTO,AJ,I) 

237* 170 CONTINUE 

270* 00 170 1*1, INODE 

271* 00 170 J*l ,3 

272* IJ-37I1-II.J 

273* 170 UOl I Jl.yl t I Jl *DUC I 1 J> 

Z 77* {•••••PRINT o{|T 

275* WRITEI6.7) III.NNn 

274* *RITEIa'7) 

277* 00 ISO l«l, [NODE 

278* J«3»l-2 

279* J|7J*I 

250» J2»J*Z 

251* IRQ WRITE (4'44) I ,UOl J) ,U0( Jl ) ,UOI J2) ,TOI I ) 

2S2* OS t G I *D5 I G2 

253* 0SIG2»0, 

257* WR I TE I 4 , 70 1 

255* DO 150 1*1, NELEMT 

254* SIGMABl jl.SIGMIBI 1 ) *0 S I G 2 B ( I 1 

257* WR1TEI4.69) ! , 05 I G ZB I I » , I S I GM A I J , I ) , J- i , 4 ) , S 1 Gm A X ( I ) , S 1 6„ AB II) 

258* 150 OS I G2"DS I G2*DS I G2B I I ) • *2 

259* ERROR*ARS( SORT IDSIG2I-S8RT IDS! Gl ) ) /SORT ( OS I G I ) 

240* WRITER*! Ill, NNN, ERROR. DS1G2,0SIGI 

2*1* lFIERROP.GT.EPSS) GO TO 3500 

242* GO TO 5000 

243* 5000 CONTINUE 

247* C*»»««CONVERGfO 

245* DO 140 }•!, NELEMT 

244* SIGMABl | >«5I GM IB! ! ) *05 I G2B | I ) 

247* 140 IF ISIGMaBI I > .GT.S lGMAX I I ) > S 1 GM A X { I > • S | G M A B I I 1 

248* (•••••NEXT TIhE INCREMENT 

249* t I I «l I I • I 

270* IF I I I I , 6* , I TER 1 STOP 

27 | • NNN«0 

272* RE AD I 5 , 2 ) t OHV I I > , I • I . NELEMT ) 

273* 17 CONTINUE 

277* WR I TE I 4 , I 0 ) IDHVI I ), I«| .NELEMT) 

275* (•••••SHIFT 

274* DO 190 i«l, INODE 

277* T2II)"T|II) 

278* Tl<II-TO<I) 

279* 0T2( I ).fiT| ( l ) 

280* DTI I I ).|>TOI II 

28 I • DO 190 j-l ,3 

282* IJ»II-1|»3.J 

283* U21 I Jl.ul 1 I J) 

287* U| I I Jl.nOI J J) 

285» DU2 I I J ) »DU i I I J ) 

284* DUI! I JI.DUoM J) 

287* 190 CONTINUF 

288* DO 200 j ■ I .NELEMT 

289» SIGMIBl i )»SIGMAR! | ) 

290* DS IGIB I t ) -DSI G2B I I ) 

291* 00 200 J • I ,4 



292* 
293 * 
299* 
295* 
296 • 
297* 
290* 
299* 
300* 
301* 
302* 
303* 
309* 
305* 
300* 
307* 
300* 
309* 
310* 

3 1 1 • 
312* 
313* 
3IM 
315* 
3U* 
317* 
3 I ® •' 
319* 
320* 


S I6HA I (’j»I )•$ I6M* ( J, f ) 

200 CONTINur 

TIME-TImE*DELT 
WRITE <6, 6) I II, TIME 
GO TO 3500 

1 FORMAT(aIS) 

2 FORM A T ( 0 F 1 o , 0 ) 

3 FORMAT urio.o ; 

H FORMAT ( /9X, • i t ME STEP -*«|5,‘ t TERAT I OR N0t*«l5/> 

5 FORMAT! f5»EI 3. SI 

6 FORM AT ( 1 H 1 , * •••T ! ME INCREMENT • , I 5 » 5 X • • T 0 T At TImE •MPC13.5/) 

66 FORMAT | »t , I5.HE1 3,5) 

7 FORMAT! NODE NO.* #5X#‘U* • I 9X #• V *, 1 9X •*»* »8X , t TEMPERATURE */ ) 

6 FORMAT! | 6 I 5 ) 

9 FORMAT ! | H ,I0X,» III • *#IStl0X#* NNN - •.tSitOX** ERROR -‘#£13.5, 
•*DSIG2**,El3.S,*0SlGt-**EI3.5) 

10 FORMAT ( ^HO, *OHV« # ,8E I 3,5 I 
|| FORMAT! J2E 10* 3 » 

07 FORMATmH ,2I6»8E|3«5) 

00 FORMAT! fMO, ‘ELEMENT I T I ELO S16MABAR SIgMAX SIGMAY 

• SIGMAZ SIGMaXY 5IGMAY7 SIGMAZX S 1 GH AH AX * ) 

09 FORMATiTH ,I3,9E!2.9> 

90 FORMAT ( »M0, ‘ELEMENT DSIGMabaR SI6MAX S I GMA Y SIGMAZ 

• S ( GMa X T SIGMAYZ SIGMAZX S! 6 HA-MAX SIGMB-I*) 

500 FORMAT!.^ IPCI3.S) | 

003 FORMAT ( ^ 2 I 31 
110 CONTINUE 
I I 9 CONTINUE 
END 


END Of COMPILATION! NO DIAGNOSTICS. 





1 * 

2 * 

3* 

9* 

5* 

0 * 

7* 

0* 

9* 

10 * 

II* 

12 * 

13* 

19* 

IS* 

16 * 

17* 

10 * 

19* 

20 * 

21 * 

22 * 

23* 

29* 

25* 

20 * 

27* 

20 * 

29* 

30* 

31* 

32* 

33* 

39* 

35* 

36* 

37* 

30* 

39* 

90* 

91* 

92* 

93* 

99* 

95* 

90* 

97* 

90* 

99* 

SO* 

51* 

52* 

S3* 

59* 

55* 

50* 

57* 


SUBROUTINE 0 * T A ( Fp , BS* , C I I 

COMMON /BLkA/NELEmT, 1N00E.NB ,NBHC,LHEAT .LSTRES.nN ,NM,NNH # NmH 
COMMON /BLKB/NOPRNT , NPaGE ,L TNES , MM ,kode »ktotal ,mTOT al .MlOAD.KLOAD 
common 2bLKC/CSK!29,29»|MC!0,0) ,hSK( 8 , 0 > # T YPEH ( ft # 8 > # TYPE I 1 8 , 8 > , 

• TYPEJ(8,8> .TYPEMCB.ei ,TYPEN(B#e) ,TYPES(B,8) 

COMMON/pL 0/S*< (20000) ,HK( 2000) , CK(2000) »HCK (2000) • lBNO( 100) • 

• IHBnoISq) 

COMMON /rLK/M A ( 30 > ,MB(30) »MC(30) ,M0(30) ,MP(30) ,Mq( 30) iMR(30» ,M5(30> 
COMMON >BLKl/A(IO)|BUO),C(IO),P(33)|O(33)|R(33y,«f(0)|H<6)iIPT 
COMMON 2BLK2/H|rR)|H2!6>,H3(8),H9!S| v HSr8)|N6(8y,H7m ( HB(8)i 

• H9(0) # H y 0 ( 8 ) ,G1 (R) ,62(0) , G 3 ( ® ) ,69(8) , G 5 ( 8 ) , G6 ! 8 y » G7 ( 8 ) (60(8) 
C0MH0N/bLK9/T0(60) i T1 (60) |T2(60) » OTO ( 60 ) • OT I ( 60 ) |0T2(60) 
C0MH0N/rLK3/NS(S0»9) , SC I 50 , 3 > , ARE A9 ( SO ) 

C0HM0N/ftLK5/IGAUS| ITER»EPSSi IGE 

C0MM0N/ftLK6/D<6,&) i E II 0 ) • EP ( I 0 ) # 5 V I ELD ( ! 0) , XNU ( Jo > # TKX I 1 0 1 • 

• SH( 10) .ALPHA <101 i DENST Y ( 1 0 I , NY ( 30 ) , ryELT , R T 
C0MM0N/pLK7/X (60) ,Y (60) |Z (60) ,MN0DE(8 ) ,MKODE( 29y 
DIMENSION FRM ) • C | ( I ) # BSK ( | ) ,TFR(SO) , JM ( 50 ) # I M ( $0 ) 


READ'S, f) IGAUS 

IF I IGAUS.EO. I ) GO TO 1000 

R£AD(5|ft02)(P( I I • Q C t) »R( I ) »!*l #33) 

READ ( 5 , ft02 ) ( A ( I I .B( I ) .CU) ,1*1 1 10) 

WR|TE(6^90Q) I P I I > .1-1.331 
«RITE!6 X 900) <Q( I ) , l-l ,33) 

RRlTE(6j900) (R( H, l-l .33) 

WR ITE (6 1 900)( A( I ) # I - I • I 0 ) 

WRITE(0 X 9OO) C B « ! ) , !•! . 1 0 1 
WRITE t 6 j900 )(C(I),I-l,|0) 

READ (5,599 ) ( ( T YPEH ( I , J I » I ■ I # 8 I , J- I , B ) 

READ (5,599) ( ( TYPE1 ( I • J) • l-l # 8 ) , J-l ,0) 

READ (Sj 599) ( I T YP£ J < I » J ) , t ■ I , 8 ) , J- I , 8 ) 

READ (5,599) ( (TYPeMI l , J ) , l ■ I , 8 ) , J- 1 , 8 ) 

READ (5,599) < ( T YP£N ( I # J ) 1 I ■ I , 8 ) # J- I , 8 ) 

READ (5,599) ( (TYPES! I 1 J) .U|, 8 ), J-1,0) 

READ! 5 i g02 ) ( (ESK( I , J ) 1 I - I # 29 ) , l , 29 ) 

WRITE (0 ,900) ( (ESM !,J) #1*1 ,29) # J- | #29) 

READ (5, 80;) ( (MSK( l , j) , I-l #8) # J-l #8) 

WRITE (6,900) ( ( H5M | , J ) #I-|, 8 ) , J* I # 8 ) 

READ (5,802) ( (HC( 1 , J) » |-l ,8) • J-I #8) 

WRITE (6,900) ( (HC( I # J ) ,1 — 1 # 8) 1 J — I , 8 ) 

599 FORMAT!*! IRE) 3.5) ) 

802 FORMAT!*! IPEI3.5) ) 

900 FORMAT(p ( IPEI 3.5) ) 

1000 CONTINUE 

I FORMAT ( | 6 I 5 ) 

RE AO 15, | I nCLEMT# INODE ,NB,N 0 HC# |PT, |TER*NXE 
READ(5,9> rt.epss.delt 

RE AD (5, 9 ) IE (I) ,EP( I ) .STfELD! I ) , XNU I I ) , TKX ( I ) #Sm( 1 ) • ALPHA I I ) # 

• DENSTTI I ) ,I-| » NK£ ) 

REA0(5, J) (Ntm • !»l .NELEMT) 

WRITE (6,1 INELEMT, I NODE # Nfi # NBHC # I PT # | TER • NKE 
WRITEl6,50p) RT,EP5S,0£LT 

WRITE (6,500) (E( I ) iEP( I ) ,SYIEL0( 1 ) # XNU ( I 1 ,TKX f Iy ,SH I I I # ALPHA! I ) # 

• OENSTYH ) , 1-1 ,NKE) 

WRITE(*”| ) INV( I ) t I-l .NELEMT) 



■tp 


58* 

500 

FoRHATibI 1 p£ 1 3 « 5 ) 1 

115* 

59* 

9 

F ORMAT 1 nF 1 Q ,0 ) 

Il6* 

60* 


NE*NELEmT 

117* 

61 • 


RE AD <sl2)IMAf|t.HBU>.MC»|l,HDIII.MP||>,HalI),HRI|>.NS«|I,l-l,NEl 

118* 

62* 


WRITEIA.2> (MM l 1 .HB! l 1 ,MC»I ) .Will .HPII) .HOI It .mRC II .MS! I I . t-l.NE) 

119* 

63* 

2 

FORHATIntS) 

!20« 

64* 


READ Is" 31 I X 1 ! ) ,T| II .21 1 > • I'M .INODE) 

121* 

65* 

3 

FORHAT (^FIO.O) 

|22« 

66* 


WRITE (A ,500) <X|D, TUI, 211). t-l . tNOOE) 

123* 

67* 


CALL COmTH 

I 24* 

68* 


NOPRNT-j 

I 25* 

69* 


NOPRNT.o 

126* 

70* 


NPAGE-0 

127* 

71* 


KOOE-B 

1289 

72* 


MH»3 

129* 

73* 


NP AGE • 1 

130* 

74* 


LINES") 

1319 

75* 


NNH-INOdE 

1 32* 

76* 


NHM-NNh.NBhC 

133* 

77* 


NN"HM" InODE 

134* 

78* 


NH-NN-Ng 

135* 

79* 


READI5, | ILSTRES.LHEAT 

136* 

80* 


WRITECAjI ILSTRES.LHEAT 

137* 

• 1* 


READIS,i)KlO«O.HLOAD 

138« 

82* 


writeu,i>kloao,hloao 

|39* 

83* 


IFlLHEAT.NE.il GO TO 2S0 

1 4 0 * 

84* 


HTOTAL"(NNh«(NNH«| ) i/z 

14|9 

85* 


CALL ZEROIhR.HTOTbL.I I 

142* 

86* 


CALL ZEROICK.HTOTAL.il 

1 4 3 • 

87* 


IFIHLOAQ.EQ.O) SO TO 5S0 

I 4 4 • 

88* 

16 

FORHATI.IFlO.O.HOM 

145* 

899 

502 

FORHAT | j I IPEI3.S. 1 10) 1 

1 46* 

90* 


RE AD 15, 7 20) 1 I HBN D 1 1 ) . TFR 1 I I . 1 • I , NBHC ) 

j 47 • 

9|* 

720 

FORHATliS.FiO.OI 

1 48 • 

92* 


WRITE I A, 502) ITFRI 1) .IHBNOt | ),|.|,NBHC) 

149* 

93* 


00 901 l-l, NBHC 

150* 

949 


KX-IHBNO< i ) 

JS | • 

95* 

901 

T01 KX ) bTFr 1 I 1 

152* 

96* 


WR I TE 1 6 j AO 1 ) 

153* 

97* 

601 

FORHA T | 7 X , . NODAL INPUT TEHPERATURE VECTOR •/) 

I 54 • 

989 


»RITE(a][A0q) 1 TO 1! 1.1*1 ,NNH) 

i55* 

999 

600 

FORMAT! | PE | 3 . 5 ) 

156* 

loo* 

550 

CONTINUr 

157* 

1019 


REWIND j 

| 58* 

102* 


REWIND q 

159* 

|03« 

40 

L-l 

160* 

1049 


NX-0 

161* 

105* 

41 

CONTINUE 

I 62* 

lot* 


H-NT(L) 

1639 

1079 


CALL HC0DEIHKOOE. HNOOE, L.HNI 

I64» 

1089 


1 F 1 IGAUJ.EO.O) GO TO 1010 

165* 

1099 


CALL PORlX ,T, 2, HNoOE. NX 1 

166* 

ilQ* 


CALL CljSTc IO.EIh) .XNUIH) | 

167* 

IH* 


CALL STjrFlO.NX.NOPRNT) 

168 * 

1129 


00 395 |-1,B 

169* 

1139 


00 395 >I,B 

170* 

I 149 


HCII.J|.HC|I.JI-ShIH)-2./0ELT 

1*7 1 • 


N 


345 HSK U , J)>HSK( I v J)«TKX(H) 

1010 CONTINUE 

WRITE(3)ESK ,NK0Dt ,MN00E 

WRITE!**) TVPEH, TYPE I , T YPEJ , TYPEM , TYPEN # TYPES 

CALL KStACK!CK,HC , 8 , 1 , HNOOE , 0 , NNH I 

CALL KSTACK<HK»HSK«8 9 ! , HNOOE . 8 , NNH I 

IF<L»6E.N£lEMTI Go to 42 

L*L* I 

NX«| 

If 1NYILI •KlE.NYU-i >) NX«0 
GO TO 4| 

42 IMNOPR^T , EO • 0) GO TO 43 
CALL PRjNTk <CK ,BSK ,NNH) 

CALL PRiNTK<HtC.BSK v NNH) 

*»3 CONTINUE 

c • • • • • surf ace HEaT BOUNDARY 
L* 1 

C***«*NS«C0UNT*Cu0CK»5CfLI)*ALPHAB»SCCL#2I»0,SClLt3l*T'' 

50 READ(S^I) (NSfLt I I ,1-1 t4) , (SCUtl > »t-I .3) 

WRITE (6,52) I NS(L* II . I-l |4) , (SCCL.I I »f-l .3) 

51 FORHAT ( <* t 5 . 6 F I 0 • 0 ) 

5 2 FORMATION .4l5.3Ei5.7l 

If (NS (l 4 1 ) .EQ.O) GO TO 60 
I I *NS ( L j I > 

I2-NSIl J 2» 

I3*NS(L 1 3) 

I4«NS(l.**> 

AREA*Ml»« 0.25*ABS( I Y ( I t l«Y| I2))#(Z( Ill-Zt 14) » 

• ♦IY( m-YU4MMZt IIMII2II 

• ♦ 17 ! ill -Z(l2»)*(X(I1)-XM4M*(ZC!|»-E(|4>| f (X(tn«X (| 2 )I 

• ♦<* U I >*X( !2> >•^(11 l-Y( !*M)MXUi |.*c Hl| 0 (Ym |*fU2>) » 

A|*AREA**(L)*SC(Li I ) *4 . /9# 

A2>AREAr(L)«SC(L. 1 >•! ./9, 

A3*AREar(L|«5C(Li | I *2* /9 * 

DO 55 I f 1 i M 
DO 55 j.l ,«* 

N!»NS(LjI ) 

NJ«NS(L,J> 

I J«NNH*uJ-NJ* (NJ-I >/2-NNH*N! 

IflNI • L T • N j ) GO TO 55 
If (!#EO*Jl HK( f J ) «HK I I J I * A I 
If I IABSI l-J) tE0*2) HK ( J J 1 *HK ( I J) ♦ A2 
If < ltNE*J*AND. IABS< I-JI*NE.2I HK ( l J)*HK t I J)*A3 
55 CONTINUr 
L*L ♦ l 
GO TO 50 

60 continue 

DO 346 i*| .MTOTAL 

346 HCK( I )*cf( ! I ♦HK t I I 
DO 544 f • t .MTOTAL 

544 HK ( I l-HC* ( I » 

If IMLOA0«EQ.O>NHH.NNH 
IF(HLOAO.EO.O) 6-Q TO 545 
CALL BOuNDiHCK 9 IHBND*NNM,NBHC) 

545 CONTINUE 

CALL PR|NTk(HCKiBSK,NMHI 
CALL HATlNv(HCK.C| ,NMH| 




CALL PRjnTk (HCK » B SK * NMH ) 

229* 

CALL BOuND(SK t !BND t NN>NB) 


00 305 j ■ 1 .MTOTAL 

230* 

CALL HATlNVtSK'Cf ,NM> 

305 

HCK I I J»?.*hCK( f ) 

231* 

CONTINUE 

250 

CONTlNgp 

232* 

RETURN 


IF (LSTRrS.EO.OJ &CTURN 

233* 

END 


KTOTAL- |NN. INN* I ) 1/2 


I ■ | 

710 REA0I5.70SI NODE, |U, IV, In 
IF lNODE. E 8 . 0 ’ 60 TO 700 
1F(IU.E0*I> I BNO ( | 1 -3-N00E-2 
IFIIU.EQ.II I • I * • 
lrlIV*EO«ll IBNOI | I-3-N00E-I 

ifiiv.eo.ii ! ■ i ♦ I 

IFI IA,EO« I ) IBNOI I I -3 -NODE 
IFIIN.Eg.il 1-1*1 
60 TO 7 l 0 
700 CONTINUE 
70S FORMAT(qlS) 

■ RITE 1*^7 1 1 IBNOI I ) , T»l ,NBI 
7 FORMAT! | 7 I S I 

CALL ZEROISE .KTOTaL, I I 
IF IKLOAQ.EO.OI 60 TO S60 
READ I 5 , j II NH 
HRITEU^II INH 
II FORMATI |01SI 

C | M I | I IS the direction of the load AT THAT joint, JhUI IS JOINT number 

S03 FORMATIjl |Pfl3.S,21SI I 
17 FORMAT! j(fIo.0.2IS) ) 

read<s i i7htfrih,ihiii,jnm>,in|,nmi 

NRITEUjSOaHTFRI I ) ,IMI ! I , JMI I > , I«i ,NH| 

DO T02 | ■ I , NH 

KT-MM-jr! | |.MM*IM( I | 

702 FR(KV)-rR(xVl*TFRtI) 

SAO CONTINUF 
RENINO 3 

17 L-l 

NX-0 

20 CONTINUE 
M.NTILI 

IF ILHEaT.NE, I I 60 TO 261 
READ I 3 | ESk.MKOOE.MNOOE 
60 TO 26* 

26 1 CALL MC(jOElMrOOE.MNOOE,L,MMl 
IFI I6AUS*E0.0) 60 TO 262 
CALL P8R(X,V,Z,MN0DE,NX) 

CALL ELASTC I 0 ,E I Ml , JNUt HI | 

CALL STlFFlO.NK.NoPRNTI 

262 CALL ASSlGN(E5K,MN0DE,NX,K0DE,L*27) 

10 CALL KST»ex(SK,ESK,8,3,MK0DE,27,NNI 

IFIL.6E.NELEHT1 Go TO 30 
L-L* I 
NX- I 

IFINTIl I .N e.NTIL-i M NX-0 
60 TO 20 

30 IFINOMRnT.EO.OI GO TO 37 
CALL PR|NTk(SK,BSK.NN) 

37 CONTINUE 


END OF COMPILATION, MS DIAGNOSTICS* 


I* 

2* 
3- 
7 • 
5- 
6 • 
7- 
B- 
7- 
10- 
i i • 
12- 
13* 
IN* 
15- 
16 * 
17* 
IB* 
|7» 
20- 
21* 
*2* 
*3* 
27- 
*S- 
26- 

27- 

28- 
27« 
30* 
SI- 

32- 

33- 

37- 

35- 

36- 
ST- 

38- 

37- 

70- 

71 - 
7*o 
73- 
77* 


SUBROUTINE DF VC TR 1 EOF , K APP A , DM , NX , NNN , I E , NL ) 

COHMON/RLKq/TOI 601 ,TI (601 ,T2IA0I , D TO I 60 I ,DTI IAO| .DT2I60I 
COMMON/rLKS/IGAUS, I TER, SPSS, IGE 

COMMON/rLK6/0(6,6|,EIIOI,EPIIO),STIELOI|OI,XNU(7o),TKXIIDI, 

• SHI 101 (ALPHA! 10) .OENSTTI 101 .NVIJOI .oELT.RT 

COMMON/g LX 7/ XI 60 1 ,71601,21601 , MNODE I B ) , MX00EI27) 
COMMON/rLKB/UoI I Bo) ,UI I 1801 ,U2( | BO I .OUoUBO) ,DU| I 1801 ,0021 ISO) 
COMMON/rLXI i/STRNl6,30) ,DSTRN|6,30) , STRNP | 6 , 30 I , ABC 
DIMENSION SI (8) ,S2(B,B,8) ,53(8,8,8,3) ,S7(B,B,3) , EDF I 8 I 
NN-NTIne) 

Bt J — EInNIo ALPHA! NNI/I I , -2 , o XNU I NN ) I 
IFINX.NE-O) 60 TO 500 
IFI IGAUS.EO.O) GO TO 1000 
DO 10 1.1,8 

CALL GauSSi IO.AA , I , | , I I 
10 S I I I ) - A A 
DO 20 1*1,8 
DO 20 J.l ,8 
DO 20 X ■ I ,8 

CALL GauSS(i 2,AA,I ,J,XI 
20 S2I I ,J,x)-AA 
DO 30 I - I ■ 8 
DO 30 IK-1,3 
DO 30 J.1,8 
DO 30 K.1,8 
IFI I X • EO • 1 I IT-16 
IFI IX*E0«2) I T» | 7 
IFIIX.E0.3l IT-18 
CALL GAySS I IT ,AA, I , J , X ) 

30 531 I , J,x , IX l-AA 
DO 70 1-1,8 
00 70 Ix*l,3 
DO 70 J.l, g 
IFUK.EO.il IT. 13 
IFIIK.E0.2I IT-17 
IFI IK.EQ.3l IT.IS 
CALL GaiiSSi IT, AA, | , J, | ) 

70 S 7 I I , J.jxl.AA 
60 TO SoO 
1000 CONTINUE 

READIS, 606)51 
REAOIS, A06IS2 
READIS, 606 ) S3 
READIS, *06)S7 




05 • 
H 6 « 
07* 
* 10 * 
09* 
50* 
SI • 
52* 
S3* 
50* 
55* 
56* 
57* 
50* 
59* 
60° 
61 * 
62* 
03* 
00 * 
05* 
06* 
07* 
00 * 
09* 
70* 
71* 
72* 
73* 
79* 
75* 
70* 
77* 
70* 
79* 
80° 
ai* 
62* 
83* 
89* 
85* 
80* 
87* 
88 * 
89* 
90* 
91* 
92* 
93® 
99 * 
95* 
90* 
97* 
98* 
99 * 

Joo° 

101 * 


500 CONTINUE 

00 90 ! * I o 8 
SUM-0.0 

SUM*SUm<.S I ( 1 J *oh 
SUM-SUNjSI II )*A 8C 
DO 100 | R« | * 8 
ira-mnooE< IR» 

DO ICO |M«| t B 
IMA-MNOnE ( JH1 

C«SH(NN»*S2M » ! R » I M ) /RT /OELT 
IF t NNN.mE* l 1 5 0 TO 50 
T0I«T I < ifRA 1 
DTOI-Tj ( I K a * 

DTI2-T21 IMA) 

T I 2 * T 2 ( | R A ) 

GO TO 60 
50 TO 1 ■TQ( j RA 1 

dtoi-tommai 

0 T I 2 • T J t IMA J 
T 1 2 »T I T i A A ) 

00 CONTINUE 

S*C« (2,*TOl*0T01-TOI*DT! 2-TI 2*DT0J I 

S«C*TM !RA)MT| I tnA)~t2<!MA)) 

IM It.Eo^O. ANO.NNN.GT. M S*C*TOI I R A I ♦ ( TO ( I HA I -tT « * M A M 
100 5UH«SUH*S 
ST J*0. 

00 1 20 lK- 1,3 
5 0*0 • 

DO 130 tK*I ,8 


102* 


IFINOPRrT.EO.OI 60 TO SOI 

l03* 


IFINX.NF*0> RETURN 

100* 


WR I TE ( 0 * 70 | ) 

105* 


WR ITE (6 jTOol 1 NNODeI 1 ) il«l ,81 

100* 


MR I T E ( 0 j 005 > 

107* 


HR ITE (6 ,400) (EOFI 1 1 , !»l «8I 

108* 

701 

FORMAT ( / 1 OX . * ELEMENT NODES */) 

109* 

700 

FORMATIrIS) 

1 ! 0* 

000 

FORMATIrI IPE13.S) 1 

in* 


WR I TE ( 6 , AO 1 ) 

I 12* 

001 

FORMAT ( /20X ■ * MATRIX S| •/! 

113* 


WRITEUjAOO) ISI I 1 1 1 l>l ,B> 

1 10* 


RRITEIA 4021 

Its* 

002 

FORM A T ( /20 X * * matrix S3 •/) 

1 16* 


WRITEI& J A0ol<(IS2II>JiKliI«I t SI,J«l > e),K*l,B) 

117* 


WRITE(0 t 0O3) 

118* 

003 

FORMATl^lOX,* matrix SJ >/) 

1 !?• 


KRITEUjAOo) 1 1 HSSI I > J>R ■ IK) 1 1°r>BI • J»l ,B) ,KM[,n j 

1 20* 


WRITE (0,000) 

1*1* 

000 

FORMAT I /20X • * MATRIX SR •/) 

|22» 


WRlTEU ( BOa) III SRI t . J. IK) , l»| ,81 .J*l .81 t lK«l,31 

?23* 

005 

F ORMAT ( / 20X » v MATRIX EOF •/) 

120* 

000 

F0RMAT(*El3.S) 

125* 

501 

CONTINUE 

120* 


RETURN 

1 2T« 


ENO 


END OF COMPILATION! NO DIAGNOSTICS* 


»TK«i»S) 


1KH«IK*( |H*| 1*3 
HE«MKOOE< i k m ) 

TEHpoRT 
00 190 J Q» I » S 
IQA-MN0pE« JO) 

I P { NNN , y E • | ) GO TO 150 
TOI*TU|QA) 

OUOI-OUi (Mk) 

0UI2»0u?«Hk) 

GO TO 1^0 
150 TO I *T0 ( j OA ) 

DU01«DUg<M*) 

0U12*0Ul (M k) 

100 CONTINUE 

DUdl*Ul fMKI 
out 2»U2(MK) 

TOI*TlC !QA> 

I r I IE«EQtO. AND.NNN.GT* 1 ) DUOI-UOTMK ) 
IF< IE.EO-0. AND.NNN.GT* I ) DU I 2*U | C HK > 

ire ie*eo*o.and.nnn*gt. i » toi-touoai 

190 TEHPcVEmP*t6i/8*0 

130 S0*SQ*TeMP*S 9< I t In. I K I • < DUO 1 -OU 1 2 ) 
120 5 ! J«S I J*8 I j*SQ 

SUN«SUN*SI J*(KAPPA)/0ELT 
90 EOF < I l*sUN 


i* Subroutine dtheta <ofb,of ,nload,dto,dti ,tp» 

7* CONHON/BLkA/NELEmT|1NOOE*N8»NBMC#I.HEAT .LSTRES v nN»NM.NNH ( NhH 

3* COMMON/rLO/SK! 20000) *HK(2000) . CK I 2000 ) , HCM 2000 > • I SND < 100) * 

9® • IHBnMSQ) 

5* DIMENSION DF( I ) »OFB( I ) |0T0( I ) t DTI ( I ) *OX(0O) aOY(0O) 

0* DIMENSION TP10O) 

7* CALL MAlMPYCCK*NNH.OTI ,DFB,0O) 

fl* «RITE(0,io6)(OFB(| ) .!■! ,NNH) 

9* 00 300 i - I . INODE 

10* 300 DFB< U.pFB ( l J ♦DP < I > 

II* CALL REoUCEtDFB, I mBND.NNH,NBHC,0O) 

12* lF(ML0An*E0*0) GO TO 10 

13* DO 30 1 • I * NNH 

10* OYIM*0'0 

15* 00 30 J«1 »NRHC 


NOPRNTaj 

NOPRNToo 



16* 


NK*IHBNO< Jl 

17* 

30 

IFd.EO.NX) 0YCI1.TPCII 

18* 


WRITECAjIOO) COYC 1 1 . 1 ■ 1 ,NNH| 

|9* 


CALL H A jMP y (HR .NN h.OY , OX, 401 

20* 


RR I TE C 4 a 1 00 ICOXCII.I*! f NNH | 

21* 


CALL REpUCE COX , IHBNO, NNH, NBHC ,401 

22* 


DO 20 I.l. NNH 

23* 

20 

DfBC I 1 .pFB ( | l-DX 1 1 1 

29* 

10 

CONTINUE 

25* 


CALL ZER0(0X.40.t| 

26* 


CALL Ma T MPY(HCIC,NMH,0FB.0X,60! 

27* 


CALL REpUCE CDT 1 . I HBND , NNH , NBMC , 40 1 

28* 


00 307 1*1, NMH 

29* 

307 

0X1 I ).DX< t >»0T ICtl 

30* 


CALL RSTORjj IDX , I HBND , NNH , NBMC , 40 • 

31* 


CALL RStOR | C OT 1 . 1 MBNO , NNH f NBHC , 40 I 

32 * 


CALL RSTOR | CDFB, 1 HBND , NNH , NBHC . 40 1 

33* 


00 304 j"l # NNH 

39* 

306 

orom.pxcti 

38* 


00 JOB 1*1 .NNH 

36* 


OTOC 1 l.pTOC 114781 I 1 

37 * 

308 

OTIC I l.pTI ( I 1 4TP I 1 1 

38* 


NOPRNT. | 

39* 


NOPRNT.p 

90* 


IFCNOPRnT a EO • 0 ) GO TO SO 

9 1 * 


■RITEI«*I03I 

92* 

103 

FORHATI/. matrix OFB ’/) 

93* 


RRITECsjlOOl OFB 

99 • 


RRITEC4.I0RI 

95* 

109 

FORMAT!/' matrix OX •/! 

96* 


RRITEC4«I 00 1 C OX 1 I I . 1*1 .NNH 1 

97* 


BRITE (4.101 I 

98* 

101 

format i • matrix oti *> 

99* 


WRITEl42>OOI oti 

6Q* 

loo 

FORMAT C JOI IPE 1 J.S 1 1 

si* 


■RITEC* 1021 

82* 

102 

FORNAY({. MATRIX OTO •/) 

63* 


*R I TE C 4 . 1 00 1 OTO 

69* 

60 

CONTINUE 

55* 


RETURN 

56* 


END 


EMO Or COMPILATION! NO 0I»6N0ST|CS. 


SUBROUTINE STRSEL (OUOiDTOi I YIELD* IN) 


2* 

3* 

q* 

5* 

6* 

7* 

8* 

9* 

10 * 

II* 

12* 

13* 

19* 

15* 

16* 

17* 

18* 

19* 

20 * 

21* 

22* 

23* 

29* 

25* 

26* 

27* 

28* 

29* 

30* 

31* 

32* 

33* 

39* 

35* 

36* 

37* 

38* 

39* 

90* 

91* 

92* 

93* 

99* 

9S* 

96* 

97* 

98* 

99* 

SO* 

51* 

52* 


REAL LAMDA.MU 
REAL jr nB 

COMMON/rLK|/»CIO>iBCIOIiCCIO>>P(JJ*>8<J3>.R<J3),»(6I>H(4>.IPT 
COMMON /BLIC2/H I C 8 1 ,H2I 8 > »MJC 81 .HR I 8) .H5 I BT ,H4 t 8 | >M7 ( B I .H8 < 8 ) . 
•H*(8>iH|OI8l>GA<8),GBt8)|6CI8>,GO<e)iGE(8),SF(8i.GH(8),GI(8) 

C0MM0N/rLK 4/DC 4 .4 » .EC 101 .EPC 101 .SYIELOC I0> I XNUt lOI , T KX I 10 J . 

. SMC 10) .ALPHA! 101 .OENSTY! 101 , NYC JO) .DELT.RT 

COMMON/rLkV/X ( 80) .YC40I .2140) , MNODE C 8 » , Mr OPE I 2R 1 
COMMON /BLK9/0SIGMAC4I , SIGNAC 4, J0> .5 I GMA I C4.30I , 

COMMON/ rL X ib/SlGM*BC 30) .S 1 6MI B C 30 ) .05 16 J B ! JO I . Os l G2B ! J0» . 

• _ S1GM**C30) 

COMMON/bLK 1 I/STRNI4.J0I .D5TRNI4.J0I .STRNP’C 4 . 30 I 4 ABC 
0 CHENS I ON LAMDACBl . MU C 8 ) » PH I CB) .SSC41 .DEPSI6I .EoRCZR) 
DIMENSION DUO! I I tOTOC 1 ) 

NN«NYI Ir) 

BIJ— EcnN|.»LPN«CNN)/C|.-2.m»NUINNI I 
CALL 2EROIEDR.2R.I I 
CALL ZErOClAMOA.S.I » 

CALL ZEROCMU.B.D 
CALL ZEROIPHI .8. I I 
CALL ZERO C OS ! GMA .4.1) 

CALL ZERO I DECS • 4 » I > 

00 10 I ■ I i 28 
NK«MKODE< 1 ) 

10 EDRC I l.pUOINK) 

DO 20 I * I .8 

L AMO A I I)»H|II)*PI|>»HHII)»PI8)«HYI!).P(|7) 

MUC C I *h7 C I ) *8 1 1 > 4HR C I I.QC9I.H7C I >•*?! 171 
20 PHI III «h I ( I ) *R ( I I. HR I I )*Rt9l*H7!!>*RH7) 

00 30 I jl .8 

1 

I 2-3»!-l 
I3*J»I 

OEPSI I I.LAmOAC I I»EDRC 1 1 I»OEPSC I I 
DEPS C 2). MU | 1 )»EDR( I2)*0EPSC2) 

OEPS C3I.PMI ( I 1 »EOR ( I 3 I aOEPSC 3 1 

DEPS I R). MU ( I I.EDRI I 1 1 *L AMO* I I 1 »EOR t I 2 1 *OEps C R ) 

DEPS (S|. PM I I I ) • EOR I I2I4MUC U»EORC I JI*OEPSCS) 

30 0EPSC4I.LAM0AC I >*EDR1 !3I»PHI I I )*EORI 1 1 1»0EPSC4) 
JC0B*A(|t.pCll4BI|l«0I|l»CIl)»RCI) 

OETJ»I./CJCOB»8.* 

DO RO 1.1.4 

RO 0EPS1 I I.OEPSC I I.DETJ 

0T«0. 

DO RS 1.1.8 
MN.MNODE'II 

R5 0T.0T»0t0I MN) /8. 

00 SO I jl .4 

DO 40 Jr I . 4 

40 DS.'GMAI I 1*051 GMA I I ) »0 f I , J I.OEPSC Jl 

IFCI.LE. 3> 05 1 GMA < I I *0 5 I G M A C I »«BI J*OT 



53* 
S*» • 
5S* 
54* 
57* 
50* 
5M 
40* 
6J • 
62* 
63* 
60 * 
65* 
66* 
67* 
66* 
69* 
70 * 
?!• 
72 * 
73* 
79 * 
75 * 
76 * 
77 * 
78 * 
79* 
8p* 
01* 
02 * 
03 * 
09* 
05* 
06 * 


50 S!GMA(|'lN)« DS|GMA(I> 

SG-<S!gmA|I , IN>*S lGMA(2.1N)*Sl6MA(3,IN)l/3. 

SS( I )«St6MA( I , 1NI-SG 

S5(2)-StGMa(2,IN>-SG 
SS< 3)*SiGMa( 3, INI-SG 
SS<9)-StGMa(9,1N) 

5S(5)-SjGMa(5#IN) 

SS<6)»StGHa(6, IN) 

Sl«l .5* f SS( I l**2*SS(2»**2*SSI3>**2*2**(5Sf 0l**24.SS(S)**2#SS(6l**2| 

• ) 

S I GMAB ( f N)-$QRTfS| 1 

If (SIGMA0I IN) • GT.SIGMAXUN) ) 1 YIELD- 1 
DO 70 1 ■ I *6 
70 STRNC I ,Jn)-OEPS( t ) 

NOPRNT.j 

NOPRNT-o 

IF( N0PRnT.EQ*0> GO TO 51 
NRITE<6«01 ) 

WRITE (6, 00) IEDR( I )• 1*1 *29) 

WR 1 TE ( 6*90 ) 

WR I TE f 6 j83 MLAMO* M ) *M0« !> #PH1C !)•!■! *B> 

WRITE(6,*0) 

WRITE(6*09) 

WRITEC0J0O) (DEPS( | ) f t*1*6| 

WR I TE ( 6 * 06 ) OETJ 

51 CONTINUE _ 

00 FORMAT! f PE I 3 • 5 ) 

S| FORMAT!* LOCAL DISPLACEMENT MATRIX •) 

03 FORMAT(s( IPEI3.5) ) 

09 FORMAT!* ELEMENT STRAIN MATRIX •) 

06 FORMAT!* JCOBIAN VALUE ■ ',!PE|3«SI 
90 FORMAT ( // ) 

RETURN 

END 


END OP COMPILATION! NO DIAGNOSTICS. 



I • 
2* 
3* 
9* 
5* 
6* 
7* 
8* 
9* 
10* 
II* 
12* 
13* 
19* 
15* 
16* 
17 * 
10* 
19* 
20* 
21* 
22* 
23* 
29* 
25* 
26* 
27* 
20* 
29* 
30* 
31* 
32* 
33* 
39* 
35* 
36* 
37* 
38* 
39* 
90* 
91* 
92* 
93* 
99* 
95* 
96* 
97* 
90* 
99* 
50* 
51* 
52* 
53* 
59* 
55* 
56* 


SUBROUTINE STRSPL ( OUO • OTO # A J » I N ) 

REAL JCnB 
REAL LAMOA,MU 

COMMON /^LK | /A I I0>.B<I0)»C(|0>.P(33I,Q(33>*R(33) # W(6)*H(4>»1PT 
COMMON /BLK2/HI 1 8 ) » H2 I B > . H3 I 8 > . H9 | 0 ) , HS I 0 I *M4 I 0 ) • H7 < 0 I . MB ( 0 ) , 

• H9 I 0 > ,H|0(B) iGA(0) *GB(0) .GC<8) *GO(0) t GE<0) *GFI8) *GH!8) *0H0> 
COMMON/nLK6/0(6*6),C(lO> . EP ( 1 0 ) * S T I ELD C 1 0 1 , XNU ( \ 0 > , TK X ( | 0 > • 

• SH( 101 .ALPHA ( 10) iOENSTT! 10) . NY ( 30 > , nEL T , RT 
C0MM0N/qLfC7/X (60) ,Y(60) .ZI60) ,MNODE(B) ,MK0DE(29) 

COMMON /BLK9/DS1GMA I 6 ) ( SI GMA ( 6 «30) iSIGHA 1 1 6*30) 
COMMON/*LKi6/SIGMaB(30) |S1GMI8(30) • OS IG1B |30) ( Ds!G 28(30! • 

• 5 I GM a X I 30 ) 

COMMON/rLKiI /STRN( 6 # 30) *0STRN(6,30) ,STRNP(6*30) t ABC 
DIMENSION laMOA(0|»MU(0),PHI(8),SSC6)»DEP5(6)#EoR( 29) 
DIMENSION 0T0( 1 ) iDUOUl .AJ(6,6) 

NN*N Y I I hi ) 

BI J--E(nN) -ALPHA <NN)/( 1 • — 2 • * X NU ( NN ) ) 

CALL ZErOIEDR.29,1) 

CALL ZCrO(lAMDA. 8, I ) 

CALL ZerO(MU»B*1) 

CALL ZErO(pm1 *8,1 ) 

CALL ZErO(DSIGMAi6. 1 ) 

CALL ZEr0(0EPS,6.1 ) 

DO 10 I ■ I *29 
NK-MKODE ( I ) 

10 EOR< 1 )*oUQ(NM 
DO 20 1-1*8 

LAHOA! I ^ • H | ( 1 ) *PI 1 ) *H9( I) *P(9)*M7( I ) *P ( l 7 | 

HU( I >*H| IV 1*0(1 )*H9( ! ) *0 ( 9 | ♦HT ( I )*Q( 17) 

20 PH! ( I )-Hl< I )*R( I ) *H9 ( I ) *R ( 9 ) *H7 ( I ) *R ( 1 7 ) 

DO 30 1-1*8 

11*3*1-2 

I 2* 3 • 1 - | 

13*3*1 

OEPS! I |. LAND A ( I ) *E DR ( I I >*D£PS( I ) 

0EPSC2).MU( I ) • EOR ( !2)*DEPS(2) 

DEPS(3)«PH| ( I ) *EOR( 13) ♦DEPS ( 3 ) 

DEPS( 9 )«MU ( I ) *EDR( I I ) *LAMDA I! ) *EDR( |21*0EPs<9) 

DEPS(5)«PHi ( l ) *EDR ( I 2 ) *MU ( I )*EOR< !3)*DEP$<S> 

30 DEPS ( 6 ) *L A MD A ( I )*EDR( I 3 ) *PH I (!) *EOR ( I I »*0EPSI6) 

JCOB«A( j)*P( » )*Bt I )*0< I )*C( l )*R( I I 
DETJ*! ./( JCOB*0. ) 

DO 90 1-1*6 

90 DEPS( ! )-DEPS( I )*OETJ 
DT*0 • 

DO 95 I — I • 8 
HN*MNOOE< l » 

95 DT*0T*Dt0(MN>/8* 

00 50 1-1*6 
DO 60 Jf I t6 

60 DSIGMA(i)«nSIGHA(!)*(D(!,J)-AJ(I l J) )*DEPS( J) 

1 F < I *LE • 3 ) DSIGMA(I)-0S1GMa(!)*BIJ*0T 
50 S!GMA( I* IN) -SIGMA | ( | , |N)*DS!GMA( I) 

SG*( SIGMA ( | , IN)*S!-GMA (2* 1N)*S!GMA| 3, IN) ) /3. 

SSU )*S|GHaM »IN)-SG 
SSC2)-StGMa(2. IN)-SG 



57 * 
SB* 
S9 * 
60* 
6 ! * 
6 2 * 
63* 
64* 

6 S • 
66 * 
67 * 
68 * 
69* 
70 * 
71 * 
72 * 
73 * 

7 4 * 
75 * 
76 * 
77 * 
70 * 
79 * 
00 * 
01 * 
02 * 
83* 
04 * 
85 • 
06* 
07* 
00 * 


SS< 3 ) -S T G H a < 3. INl-SG 
SS(4)«S?GM A f4, IN) 

SS<S)-SiGMA(S, IN) 

S S < 6 > • S T G* A ! 6 . I N ) 

D$ I G2B ( j N ) «0 • 0 
00 70 1.1.3 
J* I *3 

OS I G2B ( T N) .05 1 G2P ( I N ) ♦ < 1 . 5*SS I ! > *DS T GMA C ! )*3.0*«;S(J>*DS1gMA(J) ) 
* /sigmabmni 

70 CONTINur 

00 75 I.t .6 

75 DSTRN! I * IN) .DEPSI I I 
NOPRNT. | 

NOPRNT.p 

1 F ( NOPRkjT*EQ*0 ) GO TO 51 
*RITE(6j81 ) 

WRITEC0J0O) (EOR( T ) . !■! »24) 

*RITE<6,R0) 

WRITE(6jR0) 

WRITE(6 j r 3)1LAMDA(I),HU(I),PMI(!),I.I,8) 

WRITE(6jB4) 

WRITE(6 1 P0) (OEPSf ! 1 .1*1.6) 

WRlTE<6 a 06> OETJ 

60 FORMAT! jPE| 3.5) 

61 F ORMAT ( • LOCAL DISPLACEMENT M A TR t X •) 

03 F0RMATI3* 1 PE I 3 • 5 ) ) 

04 format ( * element strain matrix m 
86 FORMATm JCOBIAN VALUE - *,|PE13.51 
90 FORMAT!//) 

51 CONTINUE 
RETURN 
END 


end of compilation: no diagnostics. 


SUBROUTINE SKAPPA 

common ;blka/nelFmt,!nooe,n0.nbhc,lheat .lstres.nn.nm.nnh.nmh 

C0MM0N/RLK6/DI6..6) .EMO) .EP< 10) .SYIELOI I 0 ) , XNU I J 0 ) . T K X ! I 0 ) . 

SM! 10) . ALPHA ( 10) .DENSTYI I 0 ) , N Y ( 30 ) , o E L T , R T 



I* 

2 * 

3* 

4* 




2 * 

3* 

9* 

5* 

6* 

7* 

e* 

9 * 
10* 
1 1 • 
12* 
13* 
19* 
15 * 
16* 
17* 
18 * 
19* 
20* 
21* 
22* 
23* 

2 9* 
25* 
26* 
27* 
28* 
29* 
30* 
31* 
32* 
33* 

3 9* 
35* 
36* 
37* 
38* 
39* 
90* 
91* 
92* 
93* 
9 9* 
95* 
98* 
97* 
98* 
99* 
50* 
51* 
52* 
S3* 
S9 * 
55* 
50* 


SUBROUTINE ST IFF (O.NX .NOPRNT) 

COMMON )blkC/ESKI29,29) , HC ( 8 , 8 ) ,HSK ! 8 , 8 ) , T YPEH ( « , 8 ) # TYPEJ (8,81 , 

•TYPEJ(8^8liTYPEM(8.8I,TYPEN(8,8I,TYPES(8t8) 

DIMENSION 0(6,6) 

DIMENSION HSK 1 ( a , 8 ) ,HC l ( 8 , 8 > 

I F ( N X • N£ • 0 ) GO TO 350 

CALL ZEpO(ESF ,29,29) 

CALL 7EPO(HSK.8,0I 
00 200 M* l , 8 
DO 200 n-M,8 
call gauss ( i , aa.m.n, t > 

200 typehim'ni.aa 
DO 210 M*1 ,8 
DO 2 10 N*M ,8 
CALL GaijSS(2,AA,M | N,1) 

2 | 0 TYPE ! ( M , N ) « A A 
DO 220 m*1.8 
DO 220 n*M,Q 
CALL GAySS(3,AA,M,N,ll 
220 TYPEJf M,N).AA 
DO 260 M* I , 8 
DO 260 n* m »8 
TYPEHlNjMJ.TYPEHfM.N) 

TYPE! (N^Ml.TYPEI (M,NI 
260 TYPE J(N.M).TYPEJ<H,N) 

DO 230 m- I , 8 
DO 230 N-I .8 
CALL GAijSS(9,AA,N,N, 1 > 

230 TVPEM(N,N)«AA 
DO 290 n * I , 8 
DO 290 N*l ,8 

CALL GAySS(5 i AA.M ,N, 1 ) 

290 TYPES ( M ,N » «AA 
DO 250 h*I, 8 
DO 250 N - I » 8 
CALL GAySS ( 6 i A A , M ,N , 1 I 
250 TYPEN(MjN)»AA 
DO 300 J ■ I * 8 
DO 300 J-l ,8 

E$K(I.JI*D(I.1I*TYPEH(!,J>*D<9,91*TY»E|M,J)*0U.6)*TYPEJ!I,J) 

E 5K I I I ,2»*TTPEHC J f J > *0 < 9 , 9 ) *T YpEH ( J i I » 

E SK I I § J * 1 6 ) * D < | ,3)*TYP£N( J ,JI*0(6,6I*TyPEN(J» I ) 

ESKU*8jj*8)*0(2»2l*TYPEM ! , J 1 *0 ( 9 , 9 ) • T YPEH U • J | *0 ( 5 . S l *T YPE J ( ! , 
ESK! 1*8, J*l6>*Df*»3>* TYPES! I tJ>*0< 5,51 •TYPES(J,J> 
ESK!!*iA,J*|6>*0l3»3)*TYPEjM*Jl*D(5,5l*TYPEni,J)* 
•D(6,6 )«tYPEHI I ,J» 

300 CONTINUE 

DO 310 T*l,29 
00 310 >1 ,29 
310 ESK( J, | )*E5K( 1 , J> 

10 FORMAT! a t IPEI 3,5) ) 

DO 9Q0 f*l ,8 
DO 900 j*l ,6 

900 HSKU,J|«MSKC!,J»*TyPEHC!,Jl*TYPEl( ! , J ) ♦TYPE J ( I , J > 

DO 270 1*1,8 
00 270 j-l ,8 



57* 


CALL GAuSS( 1 1 ,AA , I , J, 1 ) 



58* 

270 

HC( 1 , J)«AA 



59* 


noppnt * i 



60* 


NOPRNTao 



6 I • 


I F ( NOPRnT • EQ*0 I GO TO 30 



62* 


WRITE(6j60l ) 



63* 

601 

FORMAT ( / / 2 o X , 5HT YpEH ) 



69* 


• RITE (6*600) ( « TYPE MM , J) , 

I * ! 

,8) , J ■ 1 ,8) 

65* 


WRlTE(6,6J2> 



66* 

602 

FORMAT ( £ / 2qX , 5HTYPE I ) 



67* 


• R I T E ( 6 c .6 00 I ( (TYPEl(I,J)t 

l-l 

,8) , J- 1 ,8) 

68* 


•RITE(6’603) 



69* 

603 

FORMAT! //20X.5HTYPEJ) 



70* 


• RITE! 6~60o) KTYPEJUiJl. 

1*1 

,8! , J* 1 ,8) 

71 • 


• R I TE ( 6 | 6 0 9 ) 



72* 

609 

FORMAT ( //2oX , 5HT YpEM I 



73* 


• R I T E ( 6 , 60o ) ! ( T Y P E M ( I , Jl , 

1*1 

,8 I , J*l ,8 1 

79* 


• R I T E ( 6 j 6 0 6 ) 



75* 

606 

FORMAT ( //2ox . shtypen) 



76* 


• RITE (6^8 00 )UTYPeN(|,JI, 

I ■ 1 

,8 ) , J*l ,8 1 

77* 


• RITE (6 ,605 1 



78* 

6 0S 

FORMAT {//2oX,5HTYpESJ 



79* 


W R I T E ( 6 ^ 6 Oo ) ( (TYPES! I , J) , 

I » 1 

,8! , J-l ,8) 

80* 

600 

FORMAT!*! IpEI 3,5) ) 



81 * 


WRITE! 6~ ± S99 ) 



82* 

599 

F0RMAT(,H|/,5X,29HELEMENT 

STIFFNESS MATRIX/! 

8 3* 


• R I TE <6*600) < !ES* CI.Jl.I- 

1 ,29 ) , J-l ,29 ) 

8 9* 


WR ITE ( 6,701 ) 



85* 

701 

format (j[ox, • Thermal stiffness matrix •/) 

86* 


WRITE <4.61)01 1 (MS* 1 I , J> . 1* 1 ,B 1 , J«l ,8) 

87* 


•RITE(6.27l) 



88* 

271 

FORMAT ( / 1 X , ' TRANSIENT HEAT 

conduction matrix 

89* 


WR I TE < 6 , 600 ) ! ( HC ! 1 , J) » 1*1 

.8) 

, J* 1 ,8 ) 

90* 

30 

CONTINUE 



91* 


DO 900 1*1,8 



92* 


DO 900 J* 1,8 



93* 


MSK 1 ! I , j)*MSK ! I , J) 



99* 

900 

HC 1(1 , J)«HC( 1 . J> 



95* 


return 



96* 

350 

CONTINUE 



97* 


00 910 1*1,8 



98* 


DO 910 j-l ,8 



99* 


HSK ( I , J)*HSK 1 ! 1 , J) 



ioo* 

9 10 

HC ( I , J )«MC 1 < I , J) 



10 1 * 


return 



1 02* 


end 




end of compilation: 


NO DIAGNOSTICS 



I* 
2 * 
3 • 
8 * 
5* 
6 • 
7* 
8 * 
9* 

IP* 
I I • 
12 * 
13* 
I *4 * 


SUBROUTINE PST !FF(P* *PJ, IN, !GE » 

COMMON /BLk A/NELEmT , I NOOC ,NB .NBHC .CHEAT ,LSTRES #t gN ,NH ,NNH ,NMH 

common 2blKC/ESK< 2 8,28 ) • HC I 8 . 8 ) , HSK I B , 8 ) • T yPEH I 8 » 8 I , T YPE U 8 , 8 ) , 
•TYPEJ<8,8),TTPEM<0,8) ,TYPEN( 8 , 8 > .TYPES (8.0) 

C0MM0N/pLK9/0S|GH A (6 ) . S I G M a ( A » 30 ) , S I GM a I < 6 , 30 ) 

DIMENSION PKI 28,28) ,PJ( 6 . 6 ) 

DO 300 |*l *8 
DO 300 J*! .8 

PK( 1 , Jl • PjM.I )«TYPEHM , JMPJ( I ,8)*TYPEM( ! . J I *P J < I , 8 I * T YPEN II » J > ♦ 

• P j ( 8 , | ) »TYPEM( J, I ) ♦Pj<8,0)*TYpE I(I.J)*pJ( 8 t 6 )*TYPESU,J>* 

• P J ( 6 . I MTYPCNI J.l >*PJ< 6 . 8 )*TYPES< J, l> ♦* J I 6 » 6 > *T YPE J C ! • J > 
PK ( I . ) *p j ( I ,2).*TVPEM< I » J ) +P J ( 1 * 8 ) *T YPEH I I » J) ♦PJI 1 .51 •TYPENI 1 * J ) 

• *PJ( 8,2) *TYPE! ( t *J)*PJ( 8 f 8 )*TYPEM(J. I ) ♦PJl 8 . 5 ) * T Y PES < I .J) 

• ♦Pj(6»2)*TYPES(J»I )*PJ( 6 , 8 )*TYPEN(J,I )*PJ( 6 tB)*TYPEU< I.J) 


1 5 • 
16* 
IT* 
18* 
19* 
20 * 
21 * 
22 * 
23* 
28* 
25* 
26* 
27* 
28* 
29* 

30* 
31* 
32* 
33* 
38* 
3S • 
36* 
37* 
38* 
39* 
80* 
8 |* 
82* 
83* 
88 * 
85* 
86 * 
87* 
88 * 
89* 
50* 
51* 
52* 
53* 
5 8* 
55 * 
56* 
57* 
58* 
59* 
60* 
M* 
62* 


PK < I . J*| 6 )» PJ( I .3>*TYPEN< I , JI*PJU ,5)*TYPEMM , J)*PJI I , 6 ) *TYPEH 

• <!»J>*pj(9,3>*TYPES(I.J>*Pj<8,5)*TYPEl(I . J ) *P J ( 8 * 6 > * T Y PEM | J , I >♦ 

I PJ< 6 , 3 )*TvPEJU»J>*PJl 6 »S)*TYPES(J,!)*PJ< 6 ,A)*TYPENIJ,n 

PK(1*8,J*8)«PJ(2.2»*TYPEI(I , J ) *P J ( 2 , 8 ) • T V PEM ( J , I )*PJ(2,5)*TYPES 
*( I »J)*Pjl8,2)«TYPEM< ! »J)*PJ(8,8)*TYPEM( I *J)*PJ(tj»5)*TYPEN( I . J I ♦ 

1 PJl5,2l*TYPES( J. ! >*PJ<S,8)*TYPEN( J * I ) *P J ( 5 .5 > • jYPEJI I • Jl 
PK( I*8 ,j*I 6)*PJ(2,3)*TYPES( I ,J)*PJ(2.5)*YYPEl( I , J I *P J I 2 ,& ) *TYPEM 
*(J.I>*P j( 8,3>*TYPENU . J)*PJ( 8,5)*TYPEM( I . J > *P J ( 8 * 6 > *T YPEH ( 1 , J)* 
2PJ(5,3).TYPEJ(I , J ) *PJ < 5 * 5 ) • T YPES ( J . I ) *PJ I 5 , 6 ) • T yPE N ( J , I ) 
PK(!*l^J*jA>«PJ<3i3)*TYPEjU,JI*PJl3,S)*TYPES(j,! ) *P J ( 3 . 6 ) *T YPEN 

• f j , I )*Pj'C5,3)*TYPeSC! »J)*Pj(5,5)*TYPEI ( I . J ) ♦P J ( B * & ) • T Y PEM ( J . I )♦ 

2PJ(6.3)*TYPEN(|,J)*PJ(6»5)*TYPEM(I,J)*PJI6,6)*TYPEH(1,J» 

300 CONTINUE 

00 301 1*1.28 
DO 301 j-1 ,28 

301 PK < J » I )«PK ( 1 , J) 

NOPRNT» j 
NOPRNT.q 

I F « NOPRwT.NE. 1) GO TO 500 

600 FORMAT ( ft ( 1 PE I 3 • 5 ) ) 

WRlTE<6,40ij ) 

601 FORMAT ( //2oX . 5 HTYPEH/ ) 

WRI TE (6' 600) ( ITYPEHI t .Jl . l«! .8) . J* 1 . 8 ) 

WR I TE I 6 , 805 ) 

605 F0RMAT(^/20Xi5HTYpEl/l 

WRITE (6,600) ( (TYPE I C I . Jl , 1*1 >81 , 1 , 8 ) 

WR I TE ( 6,402 ) 

602 F0RMAT(//2QX.5HTYPEJ/> 

WRITE (6*600) ((TYPEJ(!»J).I*1.8),J*|,0) 

WRITEU.603) 

603 FORMAT ( ^ /2QX .SHTYpEM/ ) 

WR ! TE ( 6 , 60o ) If TYPEM (t , J 1 , I ■ | , 0 » , J* I , 8 I 
WRITE(6,408) 

608 F0RHAT<//20X,SHTVPEN/> 

WRITE (6^600 I (I TYP£N( ! , J) ,!•! , 8 ) , J*l , 8 ) 

WR I TE ( 6 , 606 ) 

606 F0RMATC//20X.5HTYPES/1 

WRITE (6, 600) ( (TYPeS(!»J).I*I.8I.J*1.0) 

WR | T E ( 6 , 6 I 0 ) 

610 FORMAT ( /2QX , • MATRIX J •/) 

WRlTEf6^4li)((PJl|,J),l*I,6).J*I.4) 

61 I FORMAT( 6 ( 1PE13.5) I 
WRITE<6,420) 

620 FORMATI//20X . • MATRIX PK •/) 

WRI TE (6 ,*0n> < IRK U » J» * I*t »28) , J«| ,28.1 
500 continue 
return 
END 


end of compilation: no diagnostics* 



I * 
2 * 
3 * 
9 * 
5 • 
6 * 
7 • 
8* 
9* 
ID* 
I I • 
I2» 
1 3* 
I *» • 
IS* 
16* 
17* 
18* 
19 * 
20* 
21* 
22 * 
23* 

29 * 

25* 

26* 


SUBROUTINE PRNTDT(A,B,X,Y,Z,M,NJ,NP,NN) 


c PRINTS OUT JOINT FORCES A N 0 DISPLACEMENTS FOR 3-D STRUCTURES llOOOOBn 

DIMENSION A (NN) ,B ( NNI • X I I ) , Y ( | 1 ,7 U I 

N P • N P • | I 1000070 

L I NFS* 1 1 1 000080 

• R I TE I 6 * I I NP l|0000 9 0 

WR I TE < 6 2 ) 1 10001 00 

DO 15 1*1 i N J “000110 

IF (LINES -*19)9.3.3 

3 LINES*) 1 I 000 } 30. 

NP* NP* 1 * * 0001 60 

MR I TE ( 6 . I 1 NP 1I000JHO 

WRITEU.2 ) 11000)50 

9 KH»M • j 11 0°^ I 70 

KHH»KH-m* I I I 000 | 80 

S WRITE U j I6)| ,X{ fl ,Y« ! ) ,ZI1| , I A I J » . J-*HH . ** ) . < 0 ( j 1 ,J«KHH,KH) 

is lines»ltnes*i 

1 5 CONTI NUf 1 1 0003 1 0 

| FORMAT ( t M) , I ISX ,5 mPA6E »I3/1 

2 FORMAT!// 5 X » SH JO J NT , 3 X » 7 HC OORO - X , 3 X , 7 hC OORD- Y # } X , 7 HC DOR D - Z , 6 X , 11000)30 

17HF0RCE-* »SX .7HF0RCE-Y.RX .7HF0RCE-Z • 10X ,7M0ISPL-X .7X , 7H0ISPL-Y , 11000390 

27X ,7MDl < p L-Z/ > _ 1 1000350 

16 FORMAT!/ *» f M.,2* f r?.3 l IK 1 F*.» l l«tr*.3 t 3« > FlP.3.l*.riO.S.II.FI0.3.l 1000)60 

|6X.lPEM«5, |X,IPE|3,S, IX, IPF13.5) I I000)7 0 

RETURN 11000930 

EN0 110009*40 


1 * 


Subroutine pr i ntf ( a ,b ,n ) 

t 0 7 OOo 6 0 

2 * 


DIMENSION All), B ( 1 ) 

10700070 

3* 


WRITE < 6” 3) 

l07000 e O 

9 * 


DO 2 I*), N 

107000’0 

5* 


WRITERS) 1 

10700)00 

6 * 


I N»N 

10700] 1 0 

7* 


J N« I 

107001 20 

8* 


DO 1 J* 1 » 1 

10700130 

9* 


B ( J ) * A ( UN ) 

10700] 90 

1 0* 


IN* I N- l 

10700)50 

1 l • 

t 

JN* JN* I n 

10700)60 

1 2* 


WR I T E ( 6 1 5 ) CR( J) , J*l , I » 

10700| 70 

1 3* 

2 

CONTINUE 

I0700J 80 

19* 

3 

FORMAT! |H1//, 2X , 33HST 1 FFNESS MATRIX (LOWER TRIANGLE) ///> 

10700190 

15* 

9 

FORMAT ( 1 2 X , 5NR0* *,19) 

10700200 

16* 

5 

FORMAT ( | X , | P 1 OE 1 3 . 9 ) 

10700210 

1 7* 


return 

10700220 

18* 


END 

10700230 


END OF COMPILATION! NO DIAGNOSTICS. 


end of compilation: no diagnostics. 


1* subroutine elastc<d,f.xnu> 


1* SUBROUTINE PR l NTT { A ,B , X . Y , z .M ,NJ.NP ,NN 1 

2* COHHON/nLKR/TOltOlfTI (601 •T2I60>|0 TO <*0> • DTI (*(>!• 0T2f 60) 

3 * C ..... PRINTS OUT NOD»L TEMPERATURES FOR 3“D STRUCTURES *•••*•*• 

9* DIMENSION AINN) t P|NN) ,X( I » »YI | 1 ,M 1 I 

S • NP*NP.| 

6* L I NE S* | _ 

7* WRITEtfc'llNP 

8* WRITE16.21 

9* DO 1 5 1 m I »N J 

10* I F I L 1 NES-99 1 R , 3 . 3 

II* 3 LINES* I 

|2* NP*NP* ] 

1)* WRlTEU'l 1 NP 

m* WRITE I 6 ,2 » 

15* 9 CONTINuf 

16* 5 WR I TE ( 6 , 16 I I ,X ( I ) , V ( | ) ,7M > . A t | ) , B( I > .T2I | I ,TI t f I »T01 t ) 

|7* I 9 L I NE5 * L I NE 5 * 1 

18* 15 CONTINUE 

|9* l FORMAT | t M| , | |5X .5HPAGE ,!3/) 

20* 2 FORMAT!// 5 X , 5h JO | NT , 3 X » 7HCOORD-* , 3X , ThCOORD- V , jX , 7HCO0R0-Z , 

21# *9X » OF* ,ToX , *OFB» . ID* . ’ T2 * t |OX,»Tl • f I OX v 4 TO* / > , 

22* 16 FORM AT I /5X , T 9 , 2X , F9 . 3 , I X , F9 . 3 , 1 X ,F9 , 3 , 5X ,5 ! I PE I ) .5 ) ) 

23* RETURN 

29* END 


NO diagnostics. 


I 1 00n50 

I I OOOp 7 n 
1 | 000080 
1 lOQQO’O 
1 1000)00 
I 1 000 1 10 

11000)30 

I 1000)60 
1 1000)90 
11000150 


I 1OO03 1 0 


1 I 000930 
11000990 


2« 

3 * 

9 * 

5 * 

6* 

7* 

0* 

9* 

10* 

1 I • 
12 * 
13 * 
19* 
15 * 
16* 
! 7 • 
18* 
19* 
20 * 
21 • 
22 * 

2 3 * 

29 * 

25 * 

26 * 

27* 

28 * 

29 * 


D I MENS I gN D ( 6 , 6 1 

C0NST*E/M|.*XNU)*(|.-2.*XNUM 

CALL ZErO ! D .6 , 6 I 

01 I * 1 »«|.«XNU 

D! 1 , 2 ) * x NU 

D! 1 . 3 1 * X N U 

D ( 2 , I l.jNU 

D( 2|2)*nU . 1 1 

D(2.3)*xNU 

D ! 3 , | 1 *x NU 

D ( 3 » 2 1 « x NU 

0(3,31*0(1.1) 

D!9 ,9 )«o«5« ( I . • 2 « • X NU ) 

0(5,5) *r>(9, 91 
0(6,6)oo(9 t 9) 

DO 20 I* I .6 
DO 20 J* 1 .6 
20 0( I , J)*fONST*D( I 1 J) 

NOPRNT* j 

noprnt»o 

I F ( NOPRnT • EO * 0 ) GO TO 3*»R 
WR I TE ( 6^ 100) 

100 FORMAT ( / 10X , I 9HMATR l X ELASTIC/) 

WRITE(6*200)I(D( 1 ,J)»!»1 .6) , J • l , 6 ) 
200 FORMATIaI 1PEI3.5) ) 

399 CONTINUE 
RETURN 
END 


ENo OF COMPILATION: 



end of compilation: 


NO DIAGNOSTICS 


|« SUBROUTINE 2ERO<A.M,N) 

2« DIMENSION Aft) 

3* K»M*n 

DO 10 laliK 
S* I 0 A ( I » *0 . a 

RETURN 

7* END 


END OF COMPILATION! NO DIAGNOSTICS. 


1 • 


subroutine reouceif.ib 

2* 


DIMENSION F ( NN J . T B f 1 > 

3* 


I H* NB 

R* 


NH*N 

5» 

1 

I-lBl IH) 

6 • 


IFI I-NH) Z.R.R 

7 • 

2 

NH I -NH-7 

8* 


DO 3 KajtNHl 

9 • 


K 1 *K ♦ 1 

10* 


3 F»K »«F(i 1) 

1 l* 

*4 

I M ■ I H“ 1 

12* 


NH*NM- 1 

13* 


IF ( 1 H.Eo»01 GO TO S 

1«4* 


GO TO l 

IS* 

5 

CONT INuf 

16* 


RETURN 

17* 


END 


ENO OF COMPILATION: NO DIAGNOSTICS. 




I* 



SUBROUTINE RSTORIIO.IB 

2* 



DIMENSION D f NN ) , I B 1 1 ) 

3* 



NH*N-NB 

*4* 



1 H ■ 1 

5* 

1 


IalBl IH) 

6* 



IF C I .GT.NH) GO TO 7 

7* 



TOR 1 *D ( 1 ) 

8* 


2 

DIM *0.0 

9 • 

3 


J* ! ♦ 1 

10* 



IF1J.GT.NHI GO TO 5 

l 1 * 



T D R 2 * D 1 j 1 

12* 


5 

Of J)«TDP» 

13* 



TOR 1-TDR2 

1 H • 



IF ( I .GE.NH 1 GO TO 9 

15* 



! - I ♦ 1 

16* 



GO TO 3 

17* 


7 

D ( I 1 “0.0 

18* 

9 


IFI IH«GF*NB ) GO To 10 

19* 



I H * I H ♦ 1 

20* 



NH*NH* 1 

21* 



GO TO l 

22* 

1 0 


CONT INUF 

23* 



RETURN 

2R* 



end 


Io 60 ao 6 n 

10600070 

10600000 

106000*0 


10600)30 
10600! *0 


10600200 
I 06002 I G 
! 0600? 20 

10600250 
10600260 
10600270 
10600200 
106002*0 
10600300 
1 06003 t 0 


END OF COMPILATION: NO DIAGNOSTICS. 


1 o 50 no 7 n 
1PS00080 
1 05000*0 
10500100 
105001 10 
10500120 
1 0500130 

I 0500 1 SO 
10500(60 
105001 70 
1 0500100 
105001*0 
10500200 
1 05002 10 


I • 
2 * 
3* 

5 • 
6 * 
7 • 
8 * 
9* 
10 * 
I I • 
12 * 
1 3 * 

1 

15* 
16* 
17* 
18* 
1 *• 
20 * 
21 * 
22 * 

2 3* 
29 * 
25* 
26* 
27 * 


SUBROUTINE CONTH 

COMMON /8LK2/H|(8),H2(e)»H3<8)»H<H8)tH5(M ( H6(8t»H7(B)»H8(8)» 
•M9<R)|HjO(81»GI(0>iG2(8).G3(a) l GH(8|,G5(8).G6(8)«G7(8) l G0(B) 
HI ( 1 

H 1 ( 2 ) ■ I . 0 

Ml (3>*l t 0 
Hit 4 *) *-1.0 
HI < S > •-£ .0 
HI ( A ) * | .0 
HI < 7 ) • 1,0 
HI ( 8 >•-! *0 
H2 ( 1 ) ■ 7*0 
H2< 2)*-|.0 
H2 I 3 ) « J .0 
H2 1 «4 ) ■-) .0 
H2 ( 5 ) » I .0 
H 2 < 6 >*-| .0 
H2< 7 )■ 1,0 
H 2 ( B > * - | • 0 
M 3 < 1 > - |.0 
H3< 2 ) ■- 7 . 0 
H3<3)*-i .0 
H3 1 *4 ) * 7*0 
H3 (5 >*-7»0 
H 3 ( 6 ) a 7 • 0 

H 3 l 7 ) ■ 7*0 

H3 < 8 ) "-7 * 0 



26* 

HS< I >*-| *0 

85* 

GS(S)“-| ,0 

29* 

HS < 2 ) | *0 

86* 

GS ( 5 ) ■ T.O 

30* 

HS < 3 ) ■ T .0 

87* 

GS ( 6 ) ■ i*0 

31 * 

HS ( S ) * J.O 

86* 

GS ( 7 ) ■ i . 0 

32* 

H S ( 5 ) •- l *0 

89* 

GS ( 8 ) • |»0 

33* 

HS ( 6 ) “-1*0 

90* 

G5( 1 )■ 2*0 

3S* 

HS ( 7 ) “ 1*0 

91* 

G5( 2)“-|.0 

35* 

HS ( 8 ) » |*0 

92* 

G5 ( 3 ) ■ 2*0 

36* 

H5 ( 1 I ■ ]*0 

93* 

65 ( S ) •- | *0 

37* 

H5<2>“ J* 

9 S • 

65(5)“ J.O 

38* 

H5(3>“-|.0 

95* 

65 ( 6 ) ■- i .0 

3?* 

H5 < S ) • 

96* 

65(7)“ J.O 

so* 

H5 ( 5 ) *0 

97* 

65 ( 8 ) ■- J • 0 

SI* 

HS < 6> • 

98* 

66 ( 1 ) ■ |»0 

S 2 • 

H5<7>« 1*0 

_ 9 9 • 

66 ( 2 ) » J.O 

S3* 

H5 < 8 ) ■ J « 

I 00* 

G6 ( 3 ) ■- i .0 

ss* 

H7 < 1 ) »-|*0 

1 0 t • 

66 ( S )«-J.O 

S5* 

H7 < 2 I • 0 

J102* 

G6 (5 >“-2«0 

S6* 

H7(3)*-|*0 

103* 

G6(6)“-x*0 

S 7 * 

H7(S)«-| .0 

I0S* 

66 ( 7 ) “ 2*0 

S 8 • 

H7<5)« J *0 

1 05 • 

66 ( 8 ) ■ | .0 

S9 • 

H7 ( 6 ) • J. 0 

1 06 • 

67(1)* 2«0 

50* 

H7 ( 7 ) ■ |*0 

107* 

67 ( 2 ) “>2 • 0 

51* 

H7C8I» ltO 

108* 

67 ( 3 ) ■- j .0 

52* 

H I 0 ( I >»-l *Q 

i09* 

6 7 ( S ) ■ J*0 

53* 

H10I2I- 1*0 

J 10* 

G7(5)“«2*0 

5 S • 

Hiom-.i #6 

!»*• 

G 7 ( 6 ) ■ 2*0 

55* 

H 1 0 ( S ) « ! • 0 

1»2* 

G7 ( 7 ) • 2*0 

5 6* 

H I 0 ( 5 > • 1*0 

113* 

67 (8 )“-| *0 

57* 

H|0( 6)*-l *0 

IIS* 

68 ( 1 )“-J*0 

58* 

H 1 0 ( 7 ) ■ 1.0 

115* 

G8(2)“ 2*0 

59* 

H|0(8)a.l *6 

116* 

G8(3)“-|.0 

60* 

DO 10 I b 1 * 8 

.117* 

G 8 ( S ) b |.o 

61* 

H6 ( I ) “H? ( I ) 

I I 8* 

G8 ( 5 > b J ,o 

62* 

H8 < I ) “H3 I I ) 

1 19* 

G8 (6 )“-| *0 

63* 

10 M9( I >“H5< I > 

i 20* 

60 ( 7 ) • J.O 

6S* 

DO 20 lal *8 


68 ( 8 ) ■ - | «0 

65* 

20 6! < I > “ J»0 

122* 

return 

A 6 * 

G2( 1 >"-i*0 

123* 

END 

67* 

62121“ 1*0 



68* 

62 ( 3 ) ■ 1*0 



69* 

G2(S)“-j*0 

END OF 

compilation: no 

70* 

G2(5)*-| *0 



71* 

62U)“ | tO 



72* 

G2 ( 7 > ■ J.O 



73* 

62(8) “-1*0 



7S • 

63 ( 1 )“-| * 0 



75* 

63 ( 2 ) ■- J • 0 



76* 

63(3)“ |*0 



77* 

63 ( S ) * |*0 



70* 

63 ( 5 ) ■- J • 0 



79* 

63 < 6 1 | • 0 



80* 

63(7)“ J.O 



81* 

63(8)“ |*0 



82* 

GS ( I >“-2*0 



83* 

GS ( 2 ) ■- j .0 



8 S * 

6 S ( 3 > “-T *0 




DIAGNOSTICS* 





1 • 


SUBROUTINE MATKPY 1 A ,N,X , Y ,NN) 


2* 


dimension A ( 1 1 , X ( NN ) , Y 1 NN ) 


3« 


00 9 I • J , N 

1 0 3Gg 1 »0 

9 • 


N A ■ | 

10300 l 20 

5* 


N B ■ N - | 

103001 30 

6 • 

1 

Y ( I 1 ■ u . 0 


7* 


DO 9 J» | , N 

10300160 

8* 

5 

Y( ! )*Y( j }*A(NA)*X(J! 


9 • 


IF ( J-i j 2*3,3 

10300190 

10* 

2 

NA»NA*NR 

10300200 

1 1 * 


NB»NB- 1 

10300710 

12* 


GO TO 9 

10300220 

13* 

3 

N A ■ N A ♦ 1 

10300730 

1 9 • 


CONT I NUE 

I03002**0 

IS* 


return 

1O3OL2S0 

16* 


END 

10300260 


END OF COMPILATION! no diagnostics. 


1 A 


SUBROUTINE BOUND! At IB ,N,I 

2* 


dimension A ( 1 ) 1 1 B ( 1 ! 

3* 


I H ■ N 0 

9* 


NH»N 

S* 

I 

(■IB! I H ) 

6* 


J- 1 

7* 


KA-I 

8* 


IF! I *LF • 11 GO To 6 

9* 

2 

KD«( ! 2*nH-j*I l*J> /2*f I-J 

10* 

3 

K C ■ F A ♦ J 

1 I* 


A ( KA ) ■ A f *C ) 

12* 


1FIKA ,r,e* IKD-Jl) GO TO 

13* 


K A ■ K A • | 

19* 


GO TO 3 

IS* 

9 

IF ( J ,gE* (1-1)1 GO TO 

16* 


K A ■ K A * 1 

17* 


J* J* 1 

18* 


GO TO 2 

19* 

S 

IF 1 J • GF • (NH-I ) ) GO TO 

20* 


K A»K A ♦ J 

21* 

6 

K D* ( NH • | N H ♦ 1 ) ) /? 

22* 

7 

kc»ka*nh 

23* 


A ( K A ) * A ( XC ) 

21* 


IFlKC .gE* FO! 60 TO B 

25* 


K A»K A* J 

26* 


GO TO 7 

27* 

fl 

IFIIH , | E • 1) GO TO 9 

28* 


1 H« IH- l 

29* 


NH ■ NH- 1 

30* 


GO TO | 

31* 

9 

CONT INUF 

32* 


return 


I 080no6n 
10800070 
108000 8 0 
1C0000 9 0 
10800(00 
10800110 
10000120 
I 0800 | 30 
10800 1 **0 
10800 1 SO 
10900160 

108001 70 
(08001 80 
10800(90 
10800?00 
10800710 
10800220 
10900230 
108002*10 
108002S0 
1 0800260 
10800270 
1 0800280 
108002*70 
10800300 
I 08003 » 0 
10800320 
1 0900 3 30 
108003**0 
10 * 003*0 
10800360 
I 0800370, 


33* 


END 


108003*0 


ENO OF COMPILATION: NO DIAGNOSTICS. 


I • 
2 * 
3 * 
9* 

5» 
6 * 
7* 
8 • 
9* 
10 * 
1 !• 
12 * 
13 * 
!*♦• 
IS* 
16 * 
17 * 
18 * 
19* 
20 * 
21 * 
22 * 
23* 
29* 
25 * 
26* 
27* 
28* 
29* 

30* 

31 * 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE mat IN y ( A ,c ,N I 

INVERSION SUBROUTINE F()R SYMMETRICAL MATRICES 

replaces matrix by its inverse 

G 0 00 ONI.Y FOR *ELL CONDITIONED MATRfCES 
••••METHOD .. cROUTS REDUTION (CHOLESKYS METHOD! MODIFIED FOR 

symmetrical matrices 
A IS ThE MATRIX TO be inverted 

• •••C IS A oUMMY COLUMN OR R OR VECTOR OF S|ZE N 
N IS TwE SIZE OF THE MATRIX TO BE INVERTEO 

••••INPUT LO#ER O r UPPER TRIANGLE MATRIX AS A 5INGLE ARRAY 

NUMBER COLUMN (OR R0*> RISE, STARTING *|Th OIAGoNAl ELEMENT 

••••LIMIT On The SIZE OF MATRIX A IS SLIGHTLY MORE THAN 200 

a and c should be dimensioned as single arrays in the 

MA IN PROGRAM 


10200290 

10200620 

10200030 

I 02006*0 

10200QS0 
10200060 
10200070 
10200080 
1 0200090 
10200(00 
102001 10 
10200(20 
I 0200 1 30 
I 02001 HO 
1 0200 ( SO 

10200 I 60 
10200)70 
10200(60 
I 0200 | 90 
10200200 
10200210 
10200220 


••••NO PROVISION FOR TRANSFER OF CONTROL IN CASE OF SINGULAR MATRICES 10200230 

102002H0 


CALL M A T I N V ( * ,C ,N > 


1 02002*0 
10200760 


10200780 


dimension a ( 1 > • C I 1 I 

IXlI.Jl. ( j-t >«N-< ( J-l 1*f J-7> 1/2 ♦!-*»♦» 
DO 6 J • 7 , N 
DO 6 1 • J » N 


1 0200300 
1 02003 10 
10200320 
10200330 



32* 


JX*J« 1 

33* 


SUN*0 

39* 


DO S K * J 1 J X 

35* 


I XX-| X< 1 ,K > 

3 A * 


JXX- 1 X { j,K> 

37* 


KXX-IXU'K) 

38* 

5 

SUM ■ 5 UH ♦ ( A ( IXXMMJXXU/MKXX) 

39* 


I X X • 1 X f f , J 1 

90* 

6 

A 1 | X X I * A ( IXXJ-SUM 

9|* 


N 8 -N ♦ 1 

92* 


Nl*N-| 

93* 


DO IS J ■ 1 » N 1 

99* 


JXX-|X( j, J| 

95* 


Cl J 1 ** 1 ,o'A( JXX) 

96* 


J|«J*I . 

97* 


DO 1U | ■ J 1 |N 

98* 


SUM* 0 

99* 


1 1 ■ 1 • l , 

50* 


DO 9 K-J 1 1 1 

SI* 


I XX«I X( \ ,K) 

52* 

9 

sum-sum*a (I XX I *C<K ) 

53* 


I X X« I X ( j 1 I ) 

S9* 

to 

C< I)»-Sl]H/MlXX) 

55* 


I*N 

56* 


N0-NB-I 

57* 


DO 12 ji«2.NB 

50* 


1*1-1 

59* 


I XX-! X < j , I ) 

AQ* 


SUM-0 

*1* 


11*1*1 

62* 


KXX-lX(f 1 tl 1-1 

43* 


DO 11 Kali .N 

49* 


KXX-KXX+I 

45* 

i 1 

SUM-SUM^A 1 K XX | «C 1 K 1 

64* 

12 

C |1>*C 1 j l-SUM/Al 1 XX ) 

47* 


DO 13 I-JtN 

40* 


A < J X X ) *f ( 1 ) 

49* 

13 

JX X- JX X ♦ 1 

70* 

IS 

CONTINUE 

7 1 * 


JXX-! X( n ,N> 

72* 


At JXX1-T.0/AI JXX 1 

73* 


RETURN 

79* 


END 


ENO OF COHPlt*T|ON! 


NO D!»GNOST|CS. 


10200320 
10200350 
10200350 
10200320 
10200300 
102003^0 
10200100 
102001 I 0 

I 0200120 

10200130 
1 02001 20 
10200150 
10200150 
10200170 
10200180 
10200120 
10200500 
1 02005 1 0 
10200520 
10200530 
1 0200580 
10200550 
10200550 
10200570 
10200580 
10200520 
10200600 
132006)0 
10200620 
10200630 
10200620 
102006^0 
10200650 
10200670 
10200680 
10200620 
10200700 
10200710 
10200720 
10200730 
10200720 : 
10200780 ' 
10200760 


1* 


SUBROUTINE ASSI' 

2* 


DIMENSION ESK1N 

3* 


MM- 3 

9* 


JF!NX*EQ»0) 60 

S* 


DO 390 1*1,29 

4* 


DO 390 j* 1 ,29 

7* 

390 

ESK( 1 ,J»*XSKl I , 

0* 


GO TO 350 

9* 

330 

CONTINUE 

10* 


DO 30 Jm I ,N 

1 1 • 


DO 10 Is! iKODE 

12* 


I 1 - ! *K Of»E 

13* 


12*! l*KnDE 

19* 


K- I *MM 

15* 


L-K-l 

14* 


M-L- 1 

17* 


ximi-esk 1 1 . J> 

10* 


XCL) -ESK U| |JI 

19* 

10 

X I K ) -Esr I12.JI 

20* 


DO 20 1 - 1 , N 

21* 

20 

ESK < ! , J)*X ( 1 ) 

22* 

30 

CONT INUE 

23* 


00 90 1 • | ,K 

29* 


DO 50 J « 1 , K 0 0 E 

25* 


J 1 *J*KOnE 

26* 


J2* J 1 *K pOE 



27* K*J*HM 

20o L *K • I 

29* M-L-l 

30* X<M)*EsicU . J> 

31* X(L>nESK< I .Jl > 

32* 50 X(K)«»ESr < f , J2) 

33* 00 60 Jal »N 

39* 60 ESK M I J>"X( J) 

35* 90 CONTINUE 

36° 350 CONTINUE 

37* 00 320 I ■ I *29 

38* DO 320 j*l ,29 

39* 320 XSK ( I » Ji* E 5KM , J> 

90* NOPRNTaj 

(f 1 • NOPRNT 

<t 2 • IFINOPRnT.eO.OI GO TO 39R 

S3« WR ITE <»"41 l INE . IPXU ) » ! ■ ! .81 

99* 611 F0RMAT(/5X,|9MMEH B ER NU H B E R * , I 3 , 2 X , 3 HM A » . I 3 » 2 X ,3 HM B * . I 3 , 2 X , 3 M MC « 

950 •13.2X,3MH0-,l3,2X,3HHP-,I3,2X,3HMQB # t3,2X l 3HMR« f I 3 . 2 X , 3 M«S « , I 3 / ) 

96* *R I TE ( 6*9 ) 

97* 9 FORMAT! //9x .29MELEHENT STIFFNESS MATRIX/) 

90 * WR | TE ( t | 2 ) C CCS ** f C J ) t !■! ' N > •**»! t N ! 

99* 2 FORMAT (// 1 X ,0 ( 1PC I 3*5) ) 

SO* 3 99 CONTI NuF 

51* RETURN 

S2 * ENO 


end of COMPILATIONS NO DIAGNOSTICS. 


SUBROUTINE PQRIX»Y»Z»MX»NX> 

COMMON /BLX I /A | 10) *B 110) *C C 101 .PC33) *0| 33) ,R( 33| iR| 6 ) ,H! 6 »,!PT 

dimension xm.rm.nUiHxm 
I M NX .NE .0 ) GO TO 3S0 
CALL ZERO I W , 6 . I ) 

CALL ZER0(h.6 ( |> 

IEC IPT.iiE. 6 ) GO To 20 
« ( ) >*0. ITI329992 4 * 

R( 2) ■0.3607^)5730 
«M3)*0. 96791 39399 
W{9)«W( 3 ) 

W I &>■«<?) 

»(*»«*( I ) 

H ( I ) °0, *3296*5 I 92 
H{2)«0.^6l2b’3B * 5 

Ml 3)»0. 73066)91861 

M(9)--H(3) 

H(S)*-Hf 2) 

M I 6 ) ■•H< I > 



I * 

2* 

3* 

9* 

5* 

4* 

7* 

8* 

9* 

10* 

II* 

12* 

13* 

19* 

15* 

16* 

17* 

10* 

19* 


20* 


GO TO i 1 S 

21* 

20 

ifiipt.nE.9) go to 90 

22* 


HI 1 ) ■ 0 • r 6 | 1 363 1 1 5 

23* 


H ( 2 ) ■ 0 * 33990 10935 

29* 


H 1 3 ) ■•H ( 2 ) 

25* 


HI 9>«-hc 1 ) 

26* 


*1 I >*0. 3978598951 

27* 


* 1 2) *0. *52 l 95 1 598 

28* 


Wl 3>*Wf £ 1 

29* 


W< 9 ) • W ( 1 ) 

30* 


GO TO I 1 5 

31* 

90 

CONTINUF 

32* 


*i n *1 .n 

33* 


W< 2>*1 .0 

39* 


HI 1 >»0. 5773502691 

3S • 


HI 2) *-Ht 1 ) 

36* 

1 1 5 

CONTINUE 

37* 


1 | ■ H X 1 | > 

38* 


I 2*M X 1 2 1 

39* 


I 3"H X 1 3 ) 

90* 


I 9«MX 1 9 ) 

9 I • 


I 5 ■ H X 1 5 ) 

92* 


l 6 * M X 1 6 J 

93* 


I 7 -H X C 7 > 

99* 


i 8 «mx ( a > 

9 5* 


Al 1 ) **X f t 1 l*X( ! 2 ) *X I I 3 ) *X ( I91-xriS)*X( |6)*x(l7).X(|B) 

96* 


A I 2 ) * - X f f t 1-XI !2)*X! 1 3 ) * X ( I 9 ) • X I !$)■*( I 6 ) * X I I 7 ) ♦ X 1 18) 

97* 


A 1 3 ) *-X ( I 1 )»X| I 2 ) • X | I 3 I “X I I 9 ) ♦ x 1 ! 5 ) ♦ X I I 6 ) *X ( ! 7 ) *X ( |8) 

98 * 


Al 1 *)* Xf 1 1 )-X I I 2 ) ♦ X ( I 3 ) • X 1 I 9 ) ♦ X I 1 5 ) ■ X ( I 6 )* X 1 | 7 ) • X I |B) 

99* 


AIS)« XfM)*X(I2).XIl3)-Xl!9)«XM5)-XCl6)*X<I7) + Xll0) 

50* 


A 1 6 ) • X{II>-X(I2).X|t3)*X(m)aX{15)*X!|6)*X(I7)-X(lB) 

51* 


A 1 7)-A| a) 

52* 


A ( 8 ) ■ A ( 9 ) 

53* 


A 1 9 ) ■ A ( 5 ) 

59* 


A(lO)«-xMi)*XII2)-xM3)*X|!0)*XMS)*XlI6)f*X(l7)-X{18) 

SS* 


Bl 1 >*-Yf 11 ) ♦ Y 1 1 2 ) ♦ Y | I 3 ) » Y | | 9 ) -Y 1 !S)+Y(I6)+Y( I 7 ) »Y ( |8) 

56* 


8 1 2 ) *-Y f 11 )*Y( |2)*Y|I3)*YI f 9 ) -»Y < ! 5 )-V 1 j 6 ) ♦ Y 1 !7)*YI | 0 ) 

57* 


B 1 3 ) ■• Y ( I 1 ) - Y ( I 2 ) - V ( I 3 ) • Y ( J9)*Y( I 5 ) ♦ Y 1 j A ) ♦ Y 1 J 7 ) ♦ Y 1 |8) 

58* 


BIS)* Y( I 1 1 • Y C 1 2 ) ♦ V 1 ! 3 > • Y 1 | 9 ) * Y ( I 5 ) - Y I I 6 ) ♦ Y < 1 7 ) • V 1 (6) 

590 


BI5)« Yf M > ♦ Y I I2)-YM3)«Y| I9)-YM5)-Y< | 6 ) * Y ( 1 7 ) * Y ( ! 8 ) 

60* 


B 1 6 I ■ Yf I 1 ) • Y ( I 2 ) - Y ( I 3 ) ♦ Y ( J 9 ) • Y 1 I 5 ) ♦ V 1 I 6 ) ♦ Y 1 1 7 ) .V I I 8 ) 

61® 


Bl 7 1 *8 f $) 

62* 


0 ( 0 ) «B ( q ) 

63* 


B|9)«Bl5> 

69* 


BflO)»-YUl)*Y<l2>-Yt!3)*Vt!9)»rf!S)-Y(!6|*YfI7)-Y(J0) 

65* 


CM)»-Zfn)*Zll2)*Z(13l-Zf|9)-ZClS)*ZI!6)*Z(!7)-ZCI0) 

66* 


C(2)--Z( I 1 j.ZI I2)*ZII3)*ZI |9)-Z( 15)-ZU6)*ZU7)>ZC 10) 

67* 


C I 3 1 «»Z | I 1 )-Z( 1 2 ) - Z 1 J 3 ) *Z 1 I 9 ) *Z f 151*71 | 6 > *Z < 1 7 ) +Z I |B) 

68* 


CIS)* Zim-ZllZ)*ZI!3)-ZI|9)*Zlt5)-7l!0)*Z(l7)-ZI|e) 

69* 


CI5)« Z f I 1 )*ZI I 2 ) - Z 1 I 3)»Z( |9)»ZI 1 5 > *Z 1 1 6 ) *Z 1 1 7 ) «Z I |0) 

70* 


C < 6 ) * ZftM.ZII2).ZI!3)*7lt9)-Z(I5)*ZI|6)*ZII7)-ZII0) 

71* 


C f 7 > «* C f a > 

72* 


C (8 ) «C 1 4 ) 

73* 


C(9)«Cc%) 

7 9* 


Cl 1 0 1 * - 7 1 I 1 >*Zl 1 2 ) • Z 1 1 3 ) ♦ Z | ! 9 ) ♦ Z I 1 5 ) • 7 C 1 6 1 ♦ Z 1 1 7 ) • 7 ( 18) 

75* 


DO 10 1.1.9 

76* 


A It) *0. T 25* A u 1 



77* 




70* 
IT 
80* 
SI* 
82* 
83* 
80* 
8S* 
84* 
87* 
88* 
89* 
90* 
9 I * 
92* 
93* 
99* 
95* 
94* 
97* 
98* 
99* 

ioo* 
101 * 
| 02 * 
103* 
1 09* 

ios* 
io4* 
107* 
1 08 • 
109* 
110* 
ill* 
112* 
in* 
110* 
In* 

) l 4* 

in* 

118* 
in* 
120 * 
121* 
122 * 
123* 
J29* 
125* 
124* 
i27* 
I 28 * 
129* 
130* 
131* 
I 82* 
I 33* 


toco >«0. (25. Cl I ) 

PI I t«B< ?I*C*3)-CI21*BI3> 

P(2I*B(J>*CI4)-CI2I*BI4I*B(8I*C|3)-CI81*8|3I 

P(3)*B(7)*C<9|-C<2)*8<9> 

P(9).B(0l*CI3l-C(5l*B(3) 

P(5I*B|0|*C(»I-CI 81*8(41 

P(4)*B(j)*C(IO>"6(2)*B<IO)*8(8)*C(9l-C(8l*BC9) 

P(7l*B(4l*C(9)-C(5>*BI9l 

P(8l«B(0l*CI4»-C(5»*B(4l*BIIOI*C(3l-C(|OI*BI3l 

P(25l*B(5l*C(IU»"C(S)*BIIO)*B(IO>*CI9l-C(IO>*B|ul 

PI24I*B(01*C( |0)*CI8)*B(I0> 
P(27)»B|I01.C(4)-C(I0»*B(4) 

P ( 9 1 «C ( | 1*8(31*61 1 1 *C I 3 1 
P( I0I*C ( I I *B(4I-B ( I I *C ( 4 1 

P(II)*C(I>*B(9I-BII)*C(9)*C(9)*B(3I-B(0I«C(3) 

P( I 2I*B(3I .C ( 7 l-C (3 I *8(7 I 

P( I3)*C|9)*B(9)-B(0|*C(9) 

P(I0I»C (II *B(I0I-C(10»*B(|)*C (01*8(41 -B **•»•€(*» 

P(I51»C|7I«B(9I-B(7)*C(9)*C(I0>9B(3)-B(I0)*C(3> 

P( 141*0(7). 8(41-8(71*6(4) 

P(28)*C|7).B(I0I-B(7I*C(I0)*C(I01*B(4I-B(1OI*C(a) 

P(29I»C(9I*B(10I-8(9)*C(I0) 

P(30)»C(IOI.B(9)-B(IO)*C(9) 

P( 17 )*B( I l*C (2 »-B( 2)*C ( I ) 

P( I8I*B( I I • C ( 8 1 -C ( I 1*6(8) 

P( l9)*B(0l*C(2l-C(HI*BI2t 

P( 20 )»B(l»«C(S>-etl)*8<S)*BI7)*C<2) -0(71*6(2) 

P(21 l*B(71.C(5l-C(7)*B(SI 
P(221*B(0>*CI8l-C(9l*B|8t 

P(23I*B(01.C(5I"C|0I*B(5I*B( I 0 ) *C ( 2 I -C ( I 0 ) «B I 2 ) 
PI20|*B(II*C(IOI-C(II*B(IOI*B(7I*C(8I-C(7)*B(8) 

P(3I l*B(0)*C( I u I - C ( 9 1 • B ( IO)*B( 10)*C(B)-C( 1 0 ) *8 ( n ) 
P(32I*B(7I.C(1oI-C(7I*B(10I 
B(33)*B(IOI.C(5)-8(5l*C(IOI 
S(ll*»lil*f I2l-ll2l*t(3l 

8I2)9*I3>*C|8I«*(B)*C(3»*»(4)*C«2)-»I2>*C(4) 

a(3)*»(«)*C(2l-*<2)*C(9) 

8(0>*AI3>*C(SI-A(5)*C(3) 

0(5)«»(aI*C(8)-»(8I*C(4I 

S(4I»C|7)**(I0)-‘(2I*C(I0I*C(8I*»(91-«(8I*C(9I 

Q(7l»*(9l*C(SI-««SI*C(9) 

Q(8)*C(S>*A(4>-A(5>*e(A>*C(l01*A13)-A(|0)*C(3l 

8(2SI*Cf5).AlU))-A(5l*C(tO)*CIIO)*A|9)-Atlo)*Cto) 

<3(24l*C(81.»(IO>-A<8>*C(IO) 

0(27l«CrlOI*»(4)-A(IOI*C(4| 

8 ( 9 ) *A ( T * *C (. Jl *C I I I** (3 I 
9( IOI°A< I I *C ( 4 I ”C ( 1 I »» ( 4 1 

Q(II»*»(I)*C(9I-C(II*A(9>*»(0)*CI3I-C(0I*»(3I 
8( 121=4(71. C(3I-C(7)*A(3I 
8(l3l**(9l*C(9)-t(0l*A(9) 

8(10I*A|II*C1IO>-C(I>*A(1OI-»(0I*C(4I-C(9I*»I4I 

Q(I5I*A(7I»C(9I-C|7I*A(9)*»(|0>*C(3I-C(IOI*A(3) 

SI I4»*A(71.C(4I-C(71*A(4I 

8(28I»A(7)*C(I0)-C(7)*A(I0)**II0>*CI4»-CI|0I*A(*I 

8(29I»A»0>*CIIOI-C(9I*A(IO) 

Q ( 30 ) " A ( I0)*C(9)-C( ’ 0 1 • A ( 9 | 


130* 

I 3B* 
134* 
137* 
138* 
(39* 
100 * 
J01* 

1 0 2* 

I 0 3 • 

I 00* 

[95* 

I 94* 
107* 
108* 

| 99* 

I 50* 

1st* 

I S 2* 
153* 

I 50 • 
155* 

I 54* 
157* 
158* 
159* 

| 40* 

1 4 I • 

I 42* 

I 43* 
J40* 
145* 

I 44* 
I 47 • 
I 4 8 • 
149* 

no* 
i 7 1 • 

I 72* 

in* 

I 70* 
175* 
1 7 4 • 
177 * 
I 7 8 • 
[79* 
[80* 
1 8 | • 
I 62* 

I 8 3* 
180* 
I 85* 
I* 8 4 • 
I 87* 
| 88* 
I 89* 
I 90 • 


Q( I7)*A(2)*C( I » -» ( ll*C(2l 

0 ( 18 l»A( 8 I.C( I )-A( I l*C( 8 l 
0 ( I 9 )»A| 2 )*C( 9 )-»( 9 I*C( 2 ) 

0 ( 20 I*A( 2 I*C( 7 I-C( 21 *A( 7 I«a( 5 )*C(|)-C( 5 )*A(II 
8 ( 21 I*AiS 1 »C( 7 I-A( 7 I*C( 5 ) 

Q( 22 I*A( 8 I.C( 0 I-A( 0 )*C( 8 ) 

O( 23 l»C( 9 )*A(S)-A( 0 )*C( 5 l*C(lO)*A( 2 l-A(IO)*C( 2 l 
8 ( 20 I«C(II*A( 1 O)-A(II*C(IOI*C( 71 *A 18 )-A( 7 I*C( 8 ) 
0 ( 311 »C( 0 I*A( 1 OI-A( 9 )*C<IO)*C(IOI*A( 8 )-A(IOI*C(bI 

8(321 »C( 7 l*A(l(jl-A( 7 l*C(IO) 
9 ( 33 >-C(IO)*A( 5 >-A(IO)*C(S) 

R(l)*A(j>)*B( 3 )-B( 2 )*A( 3 l 

R( 2 )*A( ? l*B( 4 )-Al 4 )*B( 2 >*A(B)»B( 3 )-A( 3 l*B(BI 

R( 3 )«A( 7 )*B( 9 )-A( 9 )* 8 ( 2 ) 

R(0)*A(bI*B(3)"A(3I*B(5I 
R ( 5 )«A (0 I *B I 41 -A (4 I *B( 8 1 

R( 41 «A( ? l* 0 (IOI-BI 2 l*A(lOl*A(Bl*B( 91 -B( 81 *A(tl 

R(7)*A(g)*BI9)-A(9l*B(5> 

R(81«A(il*B(41-B(Sl*A(4)»*(IOI*B(3)-AI3l*B(lOI 
R( 25 ) • A ( 5 ) • B I ID I * A ( I0I*B(5I*AI I0l*8|9l-B( ID)*A(o> 
R(24)*Af8).B(IOl"B(Bl*A(l01 
R( 271 *aiIOI*B( 61 'B(IOI*A (41 

R(9 I «A ( 0) *8 ( I 1-* 1 I 1 *8 (3) 

R(I0I*a|41.B(I)-»(I>*B(4) 

R(|II*Ai 31 * 8 I 01 -A( 01 *B( 31 *A( 9 I*B(I 1 -A(i)*B (91 
R ( I 21 »Ai 3 ).B( 7 )-A( 7 I*B( 3 ) 

R( I 31 «A( 91 «B( 0 |-»| 91*»(91 

R( | 01 *B( I l*A( I 0 ) - A ( I l*B( 1 O 1 *B( 0 )*A( 41 »a( 0 )*B( 4 i 
R( 15 I*B( 7 I.A( 91 -A( 7 I*B( 91 .b( 101 *AI 31 -AII 01 *B (31 
R( 141 *A( 4 ).B( 7 i-A( 71 *B( 4 l 

R(2B1«B(71.a(I01-A(71*BII0I*R(I0I*A(4I-AII0)*B(a1 

R(291«b|0)*A(IO>-A(0I*B(1OI 

R( 30 I*B(I 01 *A( 9 )-A(I 01 *B( 9 » 

R( I 7 )*A( I I *BI 2 l-*( 21 *B( 1 1 
R( I 8 )»A( I 1 * 6 (BI-B|I 1 • A ( 8 I 

R(I 91 *A( 0 >*B( 2 I-B( 01 *A( 2 ! 

R( 201 «A(I)* 8 ( 51 -P(II*A( 5 )*A( 7 )*B( 21 -B( 7 I*A (21 

R( 2 I 1 «Ac 71 .B( 5 |-B| 71 *»I 51 

R( 221 »A( 01 *B( 81 -B( 01 *A (81 

R( 23 I*A( 1 OI*B( 2 I-BIIO>*»( 21 *A( 0 )*B( 5 I-B( 0 I*»IS) 

R( 201 *A(I 1 »B(IO>“B(II*A(IO)*AI 7 I*B( 81 -B( 71 *A( 8 ) 

R( 3 l)*A( 01 *B(IOI*B( 01 *A(IO)*A(lOl»B(Bl-B(|O)*A(al 

R( 3 Z 1 *A( 71 .B(I 01 -B( 7 I*AII 0 I 

R ( 3 3 I ■ A f I0)*B(S)"B( 101 • » ( 5 1 

NOPRNT.q 

noprnt.T 

IFIWOPRyT.E8.OI 60 TO 399 
»RITE(4,30|(MXI!>.I*I.BI 

WR I TE ( 4 , I 0 o> 

*R!TE(4,20O 1(A(I), 1*1. 101 
WR I TE (4,30p> 

*RI TE(4 ,200) (Bill, 1*1. 10) 

*RITE(4,0OO) 

WRITE (4,2001 (C( I 1 ,1*1 . 101 
WR I TE ( 4,4001 

WRITE (4, 700 (P(l), 1*1. 33) 

WRITEI4.BOO) 



1 91 • 


19 2* 


J 9 3 • 


1 99* 


J95* 


|96* 

30 

| 9 7 • 

100 

198* 

200 

19 9* 

300 

200* 

900 

201 • 

90 l 

202* 

902 

203* 

600 

209* 

700 

205* 

800 

206* 

900 

207* 

399 

208* 

350 

209* 


210* 



WRITE* A^ 700 ) ( Olt » , I»l *33 > 

WRITE (6*9001 

WRIT£(6^7001 (Rill , I • I • 33 ) 

WRITE<6 1 902) 

WR|TE(6,R0il(W(I),H(n,l«|,IPT> 

FORMAT ( A I 5 ) 

FORMAT! ^tOX, 8 HMATRIX-A) 

FORMAT! | 0 ( | PE 1 3.S) > 

FORMAT! / I0X,8HMATRIX-61 
F0RMAT!^I0X,8HMATr!X-C) 

FORMAT ( ? f tPEl 3.5) ) 

FORMAT | /2Qx • • GAUSS QUORATURE constants 
FORMAT! /10x,8HMATr1X-P) 

FORMAT!*! IPE13.S) ) 

FORMAT! /IOX | 8MMATR I X-0 1 
FORHAT!/IOX,8HMATrIX-RJ 

CONTINUf 
CONT INur 

RETURN 

END 


27* 

J7« J6*Mm 

28* 

J8« J7 *Hm 

29* 

MKODE I J! 1 »N ! 

30* 

N 1 a M | ♦ | 

31* 

MX ODE I J ? 1 ■ N 2 

32* 

N2*N2* 1 

33* 

MXODE I J3 ) «N3 

39* 

N3«N3* 1 

35* 

MX ODE < jq 1 *N9 

36* 

N 9 * N 9 * 1 

37* 

MX ODE 1 *15 ) »NS 

38* 

NSaNS* | 

39* 

MKODE I JA > «N6 

90* 

N A * N A ♦ 1 

9 1 * 

MKODE < J 7 ) »N7 

92* 

N 7 ■ N 7 * | 

93* 

MKODE ( JP ) »N8 

99* 

N 8 * N 8 ♦ 1 

95* 

20 CONTINUE 

96* 

return 

97* 

END 


end of compilation: NO diagnostics. 


END OF compilation: NO DIAGNOSTICS. 


I* 
2* 
3* 
9 • 
5» 
6* 
7* 

a* 

9* 

10 * 

II* 

12* 

13* 

19* 

15* 

U* 

IT* 

18* 

19* 

20 * 

21* 

22* 

23* 

29 * 

25* 

26* 


SUBROUTINE hCODE ( MKOOE , MX , L . MM I 

dimension MKODE ( I 1 i MX ( | I 

COMMON/pLK/MA ( 30) ,MB(30> . M C(30) ,MD(30) .MP130) ,Mo(30) tMRf 30) ,MSt30i 

mx u ) “Mail i 

M X ! 2 ) "MR I L ) 

MXI3)*M£<L> 

MXl9)*Mn<l) 

MX I S 1 ■HP I L I 
MX I 6 > *M0 I L ) 

MX I 7 ) ■ MR I L ) 

MX <0 )*MS<L ) 

N| -MM* I mA I L ) - I > ♦ I 
N 2* MM • ( mP ( L ) - I > ♦ > 

N 3 * M M • ( M C I L ) - I > * * 

N 9 * M M • ( mDIl ) “ I > * 1 
N5-MM* ( mP < L ) * • > ♦ 1 
NA *MM* ( mQ I L ) * I )** 

N7-MH* ( mR I L > -I > ♦ 1 

NS*MM*(mSIl)-I I* 1 
DO 20 I a I »MH 
Jl"! 

J2" J I *Mm 
J3“ J2*Mm 
J9* J3*Mm 
J5* J9 ♦ M M 
J6* JS*Nm 


I • 
2* 
3* 
9* 
5* 
A* 
7* 
8 • 
9* 
10* 
II* 
I 2* 
13* 
19* 
IS* 
16* 
17* 
18* 
19* 
20 * 
21* 
22 * 
23* 
29* 
2S* 
26* 
27* 
20* 
29* 
3 0* 


SUBROUTINE XST ACK ( S K , C , K ODE . MM , MX 0 D E , N I ,Nn ) 
01 MENS I nN MKOOE! I).SK||).C!NI,N|),KM(8> 

NX! I .J,n)«I J-I >*N-( I J- l 1*1 J-2))/2*I-J*| 

DO 60 I a I » K ODE 
JaMM*! I-| ) ♦ | 

KaMKOOE ( J > 

60 KMlD-X 

DO 20 Km I i KOOE 
DO 20 Ja 1 »MM 

no«mm«<k-i )*J 

NK-MKODF I No ) 

N A ■ N X ( Nr *NK , NN > 

DO 20 taJiMM 
NC«MM*!r-l ) ♦ ! 

SKlNA )ai;K (NA )*C(NC iND) 

20 N A ■ N A ♦ I 
KZ«0 
KX» I 

DO 90 M.2.K00E 
K X ■ K X ♦ l 
K Z * K Z ♦ I 

DO 31 K.KX.KODE 
K Y«K-KZ 

J1 »KM( Ky I 
J 2 * K M ! K ) 

I F I J2-J| 12,9.3 

2 N I ■ J | 

NK«J2-| 

GO TO 9 

3 N J ■ J 2 



31 • 


NK» J 1 - 1 

32* 

9 

DO 31 Jal.MM 

33* 


NK * N K ♦ 1 _ 

3*4* 


NA*NX ( N i , NK ,NN 1 

35* 


00 31 I.I.MH 

36* 


I F < J! ,GT • J2> GO To 50 

37* 


ND“MM* (k-Kx ) *J 

38* 


NC-MM* (*- l ) ♦! 

39* 


GO TO | n i 

9 0* 

50 

CONTI NUF 

HI* 


ND»MM* U- l ) + J 

92* 


NC"MM*lif-KX>*! 

*4 3* 

101 

CONTINUE 

*4 9 * 


SK CNA 1 * «; K ( NA 4 *C CNC ,N0» 

H 5 • 

31 

N A ° N A ♦ I 

H 6 • 

HO 

CONTINUF 

H 7 • 


RETURN 

H 8 • 


END 


END OF COMPILATION! NO DIAGNOSTICS. 


1* 

SUBROUTINE G»USS( |T,»»,N,N,LI 

2* 

common 

/BLKI/MIO),B(IOI.C(lai,Pl33I.O(331,R(33i.*(AI.H(*l,IPT 

3* 

EXTERNAL F 

9* 

A A *0 • 0 


5* 

DO 100 

T*1 . |RT 

6* 

X a H ( I ) 


7* 

DO 100 

J ■ 1 • I R T 

8* 

t»hi j) 


9* 

DO 1 00 

r*l . JRT 

10* 

Z ■ H ( K ) 


II* 

100 AA«AA*#M)*W(J)*P(K)*F<X,Y f Z,M,N,L,!T) 

12* 

return 


13* 

end 



end of compilation: no diagnostics. 




1 • 
2 * 
3* 
H* 

s* 
6 • 
7* 

e* 

9 • 
10 * 
I 1 • 
I 2 • 
I 3 • 
1 *»• 
15 * 
16 * 
I 7* 


FUNCTION F ( X . Y . 7 »M , N , L , I T > 

COMMON /PL<2 / ' W I ( 0 > » w ? ( 8 ) . H3 ( 8 ) ,H‘M8),H5(8)»H6(8i » H7 I 8 ) .Hale) i 

•H9l8l,M70«8l*6|IA>tG?l8»,G3l8>,6 l M8>t6Sl*)*66tft)tG7f8)i6*<B) 
EXTERNAL plamoa.pmu.phi .PJC09 

G«<G|lM»*G?<M»*X*G3(M»*Y*G9<M)«7*GSCN)*X*Y*G6<M)*Y*Z*G7(MI*Z*X* 

• G84M)«x*Y*Z)/8.0 

0-(GMN)*G2fN)«X*G3lN)«Y*G‘4<N)*Z*6S<N)*X*Y + G&»N)*V*Z*G7tN)*Z*X» 

• G8fN)*Y«Y*z>/8.U 

E»(GHLI*G2{L»»X*g3(LJ»Y + G‘HL)*7 + GS(L)*X*Y*G8ILI*Y*Z*G7U.)*Z*X* 

• G8(L)«y»Y«7>/8.0 

GO TO ( 10,20,30, RO. SO. 60, 70, BO, 90, 100. MO, 120, 1 30, MO. 150, 160, 

• 1 70 , 1 8Q i , ! T 

10 FaPLAMQ A ( X , Y ,Z ,M ) .PLAMDA ( X ,Y ,Z ,N>*PJC0B< * ,Y,Z1 
GO TO 3p 0 

20 f-phu(x;y,z,m)*pmui X .Y.Z.NMPJCORIX.Y.Z) 

GO TO 300 

3 0 F-PHI (X f Y # Z # M>»PHI (X t Y i Z,N)*PJCOB(X,V t ZJ 


IB* GO TO 3oO 

|9» 90 F-PLAMoi ( X , Y ,Z ,M) *PMU< X , Y ,Z ,N) *PJC0B1 X ,Y ,Z 1 

20* GO TO 3?0 

21 • 50 F*PMUtX*. T,Z l M)«PN|(X,Y,Z,N)*PJCOR(X,Y,Z> 

22* GO TO 300 

23 • 60 F»PL AM[>a f X , Y ,Z ,M) *PM 1 (x,Y,Z,N)*PJCOBfX,Y,Z» 

2H* 300 continuf 

25* F ** F / 6 N * 

26* RETURN 

27* 7Q F.PlAMda I X , Y , Z ,M ) *0. | 25 

28* RETURN 

29* 80 FaPMU! x‘ Y ,7 ,M ) *0. | 25 

30* RETURN 

3|* 90 FmPHlI t \ Y , Z ,M ) *0» I 25 

32* RETURN 

33* 100 F*G/PJCnP ( X , Y , Z ) 

39* RETURN 

35* 110 F«G*D/Pjf OR f X , Y ,7 > 

36* RETURN 

37* 1 20 F« G*0*f/PJC0B ( X , Y , 7 ) 

38* RETURN 

39* 1 30 Fa G*PLA H OA ( X , Y ,7 ,NJ *0. I 25 

90* RETURN 

9|* 190 F* G*PMu<X f Y,Z,N)«O.I25 

92* return 

93* 150 Fa G* PhM X , Y , z ,N ) *0. 1 25 

99* RETURN 

9£* 160 Fa G*0 *pL AMD A f X , Y , Z ,L I *0. 1 25 

96* RETURN 

97* 1 70 F* G*D*pMU ( x , Y , 7 , L 1 *0« l 2S 

98* RETURN 

99* 1 80 Fa G • D • PH I f X , Y , 7 , L V • 0 • 1 25 

50 • RETURN 

51* END 


END OF COMPILATION! 


NO D 1 AGNOST I C 5 



1 • 

FUNCTION PL AMOA ( x , y ,7 ,n) 

1 • 
2 * 
3 • 

2 * 

common / 8 L k I / A ( 10) » P I 10) • C ( lO) , P ( 3 3 ) ,Q(33) , R ( 3 3 > »W(A) , H ( 6 > • I P T 

q • 

3 • 

COMMON /RLK2/H1 IP) * H2 ( 8 ) , H3 I 8 > , H9 ( 8 > ,HS I 8 ) , H6 I 0 ) »H7<8) ,H8(8) , 

s* 

q* 

•H9(fl) , hT0( 8 ) ,61(8 ) .6218) ,63(8) ,69(8) , 6 S* 8 » .6618) ,67(8 > . 68 ( 8 ) 

6 * 

5* 

EI«(p(|j*P{2)*X*P(3)*Y*P(q)*Z*P(S)*X*X*P(6)*X*Y»P(7)*Y*2*PI8)*Z*X| 

7* 

6* 

£q«(p( 9 l«P( | 0 ) • X *P ( 1 t )*Y*P( I 2 ) *Z *P 1 |3)»V*Y*P( 19)*X*Y^P( IS) • Y • Z* 

8 * 

7* 

• PI l 6 > • Z • X ) 

9* 

8 * 

E 7 ■ ( p I | 7 ) ♦ p ( ! B ) *X+P ( |9)*Y4pI20»*Z*PI2) )*Z#Z*P(2?)*X*Y*P(23)*Y*Z* 

10 * 

9 • 

•P (29 ) »Z*X ) 

I 1 • 

10 * 

EI"EI*Pf25)*X*Y*Z»P(26)*X*X*Y*P(?7)*X*X*Z 

12 * 

1 1 • 

EqBEq*p(28)«X*Y»Z-.P(29)*Y*Y«X*PC30)*Y*Y*Z 

I 3* 

12 * 

E7«E7*P(3l ) »X*Y*Z-»P( 32)*7«Z*X*P(33)*7*Z*Y 

19* 

13* 

E10*E1«y*Z+E9*Z*X*E7»X*Y 

IS* 

i q • 

E 2 *F I • Y 

16* 

15* 

E 3 “ E i • Z 

17* 

16* 

E5*E9*Z 

18* 

1 7 • 

E6*E9°X 

19* 

18* 

E 9 *E 7 ♦ X 

20 * 

19* 

E9»E 7 • Y 

21* 

20* 

PLAMDA»fMj(N)cEI*H2(N)»E?*H3(N)*F3*HqiN)*Eq*HS(N)*E5* 

22* 

21* 

*H6lN)«E6*H7(N)*E7*H8(N)*ER>H9(N)*E9) 

23* 

22* 

PLAMOA.pLAMOA4.H!0(N)*E|0 


23* 

return 



function pmu ( x,y.z»^i 

COMMON / B L k I / A < 10> *BC tO> «C| IO> .P(33)»Q(33} , R ( 3 3 > »*(6) . H I 6 ) iIPT 
COMMON ,BLK2/H1(8)»H2m,H3<8),H«4<8>,H5<P)»H6(Bl.H7(e),H8<Q»» 

• H9le).H|Ot8).6|(fi).GZ(el,G3m.69(8),G&((O,GA(8),67(8)t60<8> 

F I MQ( I ) *Q( ?) •X + Q ( 3) *Y*Q( q ) • Z *Q C S ) • X * X ♦© U ) •X*Y#Q(7 >*Y*Z*Q ( 8) *Z*X ) 
f q ■ ( 0 ( 9 ^ ♦ q ( | o ) * X ♦ Q I 1 1 I • Y ♦ 0 ( 1 2 ) * 7 *Q I |3)«T«Y*Q( |q)*X*Y*Q( |5)*Y*Z* 

• Q I i 6 ) • Z* X ) 

F7»(Qt|7>*QIIBMX*<3(|9)*Y*Qf20>*7*«(2n*7*Z*Q<22)*X*Y*Q(23)*Y*2> 

*0 ( 24 ) *Z*X ) 

Pl«F|*Qf25)*X*Y*Z*Q(26)«X*X*Y*Q(27)*X»X*Z 

FP»FM*0(28)®X*Y*Z4Q(29)*Y*Y*X40(30)*Y*Y«Z 

F7*F7*Q(3I)«X«Y*Z*QI32)*7*Z*X*Q(33>*7«Z*Y 

F|0"FI*y*Z*F9*Z*X*F7*X»Y 

F2»F I • Y 
F3»F1*2 
F5“F M * 2 
F 6 ■ F M * X 
F 8 * F 7 • X 
F 9 * F 7 * Y 

PMU»<M| ( N)*r|*H2lN)*F2*M3(N>*F3*M9IN)*FH»HS(N)*p5*HA(N)*F6* 

•M7lN)*F7*H8(N)*F8*H9(N)*F9) 

PMU»PMU4Hlft(N) *F IQ 

RETURN 


29* 


END 


29* 


END 


end of compilation: 


no diagnostics. 


end of compilation: 


NO DIAGNOSTICS 



I * 
2 * 

3 • 

*4 • 

5 >* 

6 • 
7* 
8 * 
9* 
10 * 
1 1 • 
I 2 * 
13* 
19 * 
IS* 
1 A* 
I 7 • 
18 * 
19 * 
20 * 
2 ! • 
22 * 
23* 


FUNCTION PMt(X,Y»2*N» 

COMMON / ALk 1 / A ( 101 i B ( 10 1 »C( 10) »P(33)»0(3 3) , R < 33 > i*(6) ♦ H { 6 ) ,IPT 
COMMON 2 bLK ?/n I IP) » H 2 I 8 ) » H3 » 8 ) » M*4 ( 8 ) » HS I P ) ,H6I8 j »H7l8) ,h8(8) i 
•H 9lfi),H|0IP).6AlP>*GB(8),GC<8>tGni8> # 6ElP».6F(8l.GH(8),6l<8) 

G1*IR<1»*R«2)*X*P(3>*V*P<9)*Z*R(S)*X*X*R(8)*X*Y*R(7)*Y*2*R<8)»Z*X) 

GM«IRt9i*R(!0)*X*RUI)*Y*RM2)*7*Rll31*Y*V*RllM)*X*Y*Hl|5»*Y*Z* 

• R< I 6 ) • Z • X ) 

G7«{R<l7)*pf|B)*X*R(l9)*Y*R(20)*Z+R(?I)*2*Z*RI2?)*X*Y»RC23)*Y*Z* 

*R(29|*Z*X) 

GI*G)*R|25l*X«Y*Z*R(26)*X*X*T*RI27)*X«X*Z 

GM«G8»R(28)«X»Y«74RI29)*Y*Y»X*RI30)*Y*Y*7 

G7-G7^R(3I)*X«Y»7*R<32)*Z*Z*X*RI33)*7*7*7 

G! 0 *GI*v*Z^g 9 *Z*X*G 7 *X*Y 

G2»Gl *Y 

G 3 * G 1 *Z 

GS»G*4*Z 

G6 • G ** • X 

G 8 * G 7 • X 

G 9 ■ G 7 • Y 

PHl«|H||N).Gl*H2lN)*G2»H3(N)*G3*MM(N)«G8*MSIN)«f,S*H8(N)*G6* 

•H7(NI*G7*H8IN)*G8*H9(N)*G9) 

PHl»PMl4.M|o(N)*GtO 

return 


2*4* 


END 


end of compilation: no diagnostics. 


I* 
2 * 
3* 
*»• 
5 * 
A* 
7* 
8* 
9* 
10 * 
I I • 
I 2 * 
13* 
19* 
IS* 
I A* 
17 * 
18* 
19* 


FUNCTION PjCOBIX.YiZ) 

COMMON /BLK I / a < t 0 ) i 8 ( 10) ,C( 10) , P I 3 3 ) ,0(33) , R < 33 ) .» ( A > ,H< A ) ,IPT 

COMMON >BLK 2 /H|fP),H 2 ( 8 I.H 3 l 8 ),H* 4 | 8 ),HS< 8 ),HA( 8 >,H 7 C 8 |,M 8 < 8 ), 
•H9l8»,HToi8)»GA(P)#GBI8>*GC<R>*GD|8) # GE<0».GF(8>.GMCe>,GKfi) 
E|«IP(|»*P(2)*X*P{31*Y*Pf*M*Z*PCSl*X*X*PtA)*X*Y*P(7)*T*2*P<0)*Z*X) 
E|*El*P( 2 S)»X«V* 7 *P( 2 A)*X*x*Y*P( 27 )*X«x *2 

F|»njC|)*Q42)*X*Q(3)*Y*Q(‘4l»Z*OlS)*X*X*Q<A)*X*Y 4 Q(7>*Y*Z*OI8)*Z*Xl 

F|«FI* 0 ( 2 S)*X*Y* 7 *Q{ 2 A)*X«X*Y*QI 27 )*X*X *2 

G!»CRI|)*RI2)*X*R(3)*Y*RC*4)»Z4R(5)*X*X*RIA)*X*Y*RI7)*Y*Z*RI0)*2*X) 
GI*GI*Rf25)*X»Y*Z*Rl2A)*X*X*Y*R(27l«X*x*2 
0 1 ■ A ( I ) • E | ♦ A ( 1 ) • F | ♦ C I I )*GI 
D2*A(9).E|*B(‘4)* r I*C(‘4)*Gt 

03»A(7)*EI*B(7)*F|^CI7>*GI / 

0*4*A( l0)*^l*8< 10**F|*C|10)*GI 

X JC0B»D » *D2* y *03*Z 

X JCOB-X jC0b*09 • Y*Z 

PJCOBM./XJCOS 

RETURN 

end 


NO 


end of compilation: 


D I AGNOST I CS 



